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Abstract 

C3 We calculate the most divergent contribution to non-degenerate sectorof the self-energy 

^ (or "melonic") graph in the context of the Lorentzian EPRL-FK Spin Foam model of Quan- 

tum Gravity. We find that such a contribution is logarithmically divergent in the cut-off over 
the S?7(2)-representation spins when one chooses the face amplitude guaranteeing the face- 
splitting invariance of the foam. We also find that the dependence on the boundary data is 
Q different from that of the bare propagator. This fact has its origin in the non-commutativity 

of the EPRL-FK Y^-map vnth the projector onto SL{2, C)-invariant states. In the course of the 
paper, we discuss in detail the approximations used during the calculations, its geometrical 
I interpretation as well as the physical consequences of our result. 
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1 Introduction 

The Spin Foam program is an attempt to quantize General Relativity in a background in- 
dependent and Lorentz covariant way [ 1] . Its strategy is aimed to reproduce in three plus one 
dimensions the success of the Ponzano-Regge model in providing a viable quantum theory 
of three dimensional (Euclidean) quantum gravity [2]. The result is the EPRL-FK model [3],^ 
which is an SL(2, C)- (i.e. Lorentz-) covariant [6] state sum model, defined on two-complexes,^ 
having (projected) SU{2) spin networks as boundaries. These spin networks span the same 
Hilbert space as those obtained in canonical Loop Quantum Gravity. 

Independence from the cellular complex^ (and background independence) is recovered in 
a second moment either via a refining of the complex itself or via a summation over all ad- 
mitted complexes (e.g. over all the triangulations of a given manifold). The cellular-complex 
independence of the theory is connected to its continuum limit^ which is intrinsically non- 
perturbative, and as such it is likely to require some kind of renormalization procedure; in 
turn, any renormalization procedure needs a thorough understanding of the divergences of 
the theory. 

This work aims to perform a first step in this direction for the EPRL-FK model, through the 
calculation of the most diverging part of the "self-energy" (or "melonic") cellular complex at 
fixed boundary data (see Fig. 1). This complex is dual to two four-simplices with four tetrahe- 
dral faces identified among them. The result of the calculation is that the self- energy two com- 
plex is found to diverge logarithmically in the cut-off on the SU[2) representations, at least 
for the specific face amplitude which preserves the face-splitting property of the spin foam 
ifiUf) = + !)• This result is in accordance with what was already known for the Euclidean 
version of the theory^ [9, 8] . However, here one more step is performed: the dependence of the 
amplitude from the boundary data is worked out in detail, and is found to be different from 
that of two tetrahedra directly identified among them without the mediation of any dynamical 
process. Finally, it is interesting to notice that the two spacetime parities of the virtual states 
which mediate the process play a crucial role in generating the leading divergences. 

A QFT approach to the Spin Foam program exists and is referred to as Group Field Theory 
(GET) [10] . A GET is a particular type of a QFT, whose fields live on multiple copies of a group. 
It can be shown that any spin foam can arise as a particular Feynman diagram of one such 

1 However, see [4] and references therein for a recent review of the spin foam approach to quantum gravity and 
more details about the EPRL-FK model, and also [5] for a an interesting comparison between different approaches to 
its construction. 

2 Sometimes a lower degree of generality is preferred, and the admitted cellular complexes are reduced to those 
dual to triangulated (pseudo-)manifolds; see for example [7]. 

^In the three dimensional context of the Ponzano-Regge model, this is not an issue. Indeed, modulo regularization 
difficulties, the theory is triangulation independent. This reflects the fact that three dimensional gravity is a topologi- 
cal theory, i.e. it has no local degree of freedom. 

''Remark that the continuum limit is not the same as its semi-classical limit. See scheme at p.21 of [1]. 

^Remark thatin [8], another face weight is used (i.e. f] — {2jJ +l){2jj + 1) — 17^—117^+0(77)) and the accordance 
with this paper is modulo a redefinition of the face weight. Even this looks trivial, it is important to realise that in 
the Group Field Theory model construction the face and edge weights are strictly related and cannot be modified 
independently as for Spin Foams. In the case of reference [8], in particular, the face weight is dictated by the choice of 
a "local" (trivial) vertex and a "dynamical" (non-trivial) propagator. 
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Figure 1: The vertex and edge structure of the self-energy cellular complex considered in this paper 
(see Fig. 2 for its face structure). 



GFT,^ which automatically reproduces its amplitude and at the same time provides a canon- 
ical prescription for their summation.'' In a specific class of GFTs (the so called "coloured" 
GFTs and other models derived from them; see [14, 15], as well as [16] and references therein), 
the self- energy graph studied here plays a fundamental role: in simple enough cases it was 
explicitly shown to drive the non-perturbative behaviour of such models {l/N expansion and 
critical behaviour [17, 18], renormalization flow [19, 20, 21]). This result opened the way to an 
extension of matrix models and matrix field theories to higher dimensions. In this context, the 
GET self-energy graph was first called melon graph, a terminology we shall adopt throughout 
the paper. 

Spin Foam divergences, at least in the simplest models, can be given two different inter- 
pretations, which are essentially dual to one another. The first is mainly geometrical, and as- 
sociates the divergences to unconstrained sub-simplices which are free to become arbitrarily 
large, as if some vertices of the triangulation were let free to "escape to infinity" {spikes). The 
second is in terms of unconstrained hypersurfaces within the dual complex {bubbles), which 
play a role similar to Feynman diagram loops. 

In the Ponzano-Regge model, spike divergences are known to be related to the shift sym- 
metry of BF theories and are expected to be associated to the background independence of 
four dimensional gravity [22, 23]. Still in the context of the Ponzano-Regge model, in [24] it 
was suggested that spike divergences are also intimately connected with the implicit sum over 
parities that is present in the theory.*' Here we give some arguments which suggest that the 
same might happen in the EPRL-FK spin foam model. 

The ERPL-FK model can also be extended to the case of General Relativity with cosmo- 
logical constant in a non-trivial way, both in its Euclidean [25, 26, 27] and Lorentzian [28, 27] 
versions. Such an extension uses the ^^-deformed Lorentz group, with the q'-deformation pa- 
rameter related to the cosmological constant, and turns out to be (perturbatively) finite. The 
existence of a finite model does not mean that the issue of large radiative corrections can be ig- 

^The GFT formulations of the Euclidean and Lorentzian EPRL-FK models can be found in [8] and [11], respectively. 

^Furthermore, this summation naturally implements not only a "sum over geometries", but a "sum over topologies" 
as well. For more details about GFTs, see [12, 13] and references therein. 

^The presence of a sum over parities can be tracked back to the fact that the Ponzano-Regge model, as any other 
spin foam model, stems from the Einstein-Cartan action (rather than from the Einstein-Hilbert action) which is parity- 
change sensitive. 
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nored: it may still happen that some higher order graphs have large amplitudes and therefore 
drive a renormalization flow, possibly even through phase transitions. Qualitatively, the scale 
which imposes the infrared^ finiteness of the theory is given by the cosmological constant, 
which is of the order of the radius of the Universe; therefore, at our - or smaller - length scales, 
it can be considered as infinite for practical purposes (but see comments in the conclusions). 

In this paper, in order to study the simplest EPRL-FK divergence, we introduce a cut-off A 
to the SU{2) representations j. The physical meaning of such a cut-off is that of imposing a 
maximal value for the area operator, which can be thought as the introduction of a finite size 
for the Universe itself. A bound to the area operator is typical of the ^-deformed version of 
the EPRL-FK model. Therefore the introduction of such a cut-off can be hoped to be a simple 
implementation of the main feature of the (y-deformed EPRL-FK model within the much more 
manageable non-deformed version. At the light of this (qualitative) correspondence, the cal- 
culation of this paper can be also given a more physical, though possibly naive, interpretation 
in which the cut-off A is a physical quantity and corresponds - at least in order of magnitude - 
to the cosmological constant Acc expressed in Planck units of area: A Acc/lp ~ lO^^". 

The goal of this work is to calculate the most divergent contribution to the self-energy of 
the EPRL-FK spin foam model coming from its geometrically non-degenerate sector (the ter- 
minology will be clarified later). In order to do this, various types of approximations are intro- 
duced which will be explained in detail. The main tool used is an appropriate application of 
the large spin analysis, first introduced in the context of the EPRL-FK model in [30, 31], for a 
single vertex, and in [32, 33], for a general two simplex. Thanks to the geometric content of the 
equations involved in the analysis, it will be possible to interpret the most diverging contribu- 
tion to the self-energy as stemming from the constructive contribution of two parity-related 
"chunks" of spacetime (i.e. of a quantum of spacetime with a quantum of anf i-spacetime, in 
the terminology of [34]). 

The analysis of the degenerate sector is left for future studies. However, it must be remarked 
at this point that hints coming from the Euclidean version of the model seem to say that this 
sector could be the one dominating the melon graph (see the Appendix C of [8]). An intuitive 
picture of why this may happen, expressed in the language of this paper, shall be given in Ap- 
pendix C. 

The paper starts with a brief introduction of the notation in Sect. 2. Then, in order to help 
the reader through the multitude of technical details, as well as an argument of plausibility of 
the approximation schemes, the calculation is first performed in the context of SU{2) BF the- 
ory, where the result is known and can be double-checked. This is the subject of Sect. 3. The 
definition of the EPRL-FK model on an arbitrary two complex, and the preliminary study of 
the EPRL-FK melon graph symmetries are performed in Sects. 4 and 5, respectively. Section 6 
introduces the stationary phase approximation for the melon graph amplitude and, conse- 
quently, the large spin limit of the EPRL-FK model. The equations which define the stationary 
points are then solved in Sect. 7. This will be first done in complete generality from an alge- 
braic, rather abstract, point of view, and later given a much more transparent geometric inter- 

^Here, the term "infrared" must be understood as relating to large physical distances; analogously, an "ultraviolet" 
cut-off, in the sense of a short distance cut-off, is naturally present in any spin foam models, via the existence of the 
area gap [29] . It must however be kept in mind, that the roles of the words "infrared" and "ultraviolet" are interchanged 
when referring to the momentum space of the gauge group. 
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pretation (which however is too weak to capture all the details of the algebraic solution). This 
analysis shall point out the existence of two non-degenerate sectors of solutions, related to an 
Euclidean and a Lorentzian geometry respectively, as well as a degenerate one (which deserve 
a thorough analysis on their own). The role of the two-complex face orientation is discussed 
in Sect. 8, which is crucial to correctly evaluate the melon amplitude at the stationary point. 
The radiative corrections related to the three sectors mentioned above are evaluated in Sect. 9, 
and the total graph amplitude is presented in Sect. 10. In Sect. 11, the approximations and 
hypothesis used during the whole calculations are briefly reviewed. The spike interpretation, 
as much as its limitations, is discussed in Sect. 12. Finally, we conclude in Sect. 13. 

2 Spinors, Vectors and Bivectors 

In this section we briefly introduce some notations and technical tools that will be exten- 
sively used in the following. 

Written explicitly in components, a spinor reads: 

\w) = w=^^^^eC^ and {w\=^^^^ =[w°, w^) , (2.1) 

where the bar stands for complex conjugation. In the first of these equalities we introduced 
the simplified notation for "ket" spinors. The hermitian scalar product between spinors, is 

{w\z) = w° z'^ + w'^ z'^ , (2.2) 

hence the squared norm WwW^ := {w\w} >0. 

We define the antilinear operator /: 



/ : ^ , / 



(2.3) 



It is easy to verify that {w\J\w) = Ofor all w e C^, and that that Jg J-^ = g^'^ for all g e SL[2, C). 



To a spinor one can associate a four-vector l{w): 



f 



w = 



lIwY := {w\a'\w) = 



\w\ 



22{w°w^) 
|u;0|2-|u;i|2 



, where a' = 



I 

cr' 



(2.4) 



and cr' are the Pauli matrices.^" Equivalently in a notation more adapted to normalized spinors: 



w = we 



e^"/* cos e 
e-''/'sin0 



L{ivy = \\w 



i°With the usual conventions: 



cr 





















and cr^ ■ 



f 1 A 

sin20 cos2!/i 
sin20 sin2?/^ 
\^ cos20 J 



1 

-1 



(2.5) 
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From this expression, it is readily verified tiiat u; ^ is light-like. Remark that the four-vector l[w) 
looses the information about the global phase of the spinor. The spatial part of l{uj) shall be 
noted iu, and its components ivK 



The spinor Jw is, associated to the four- vector 



I _ 



f iii^ip ^ 
-(|m;0|2-|m;1|2) ; 



(2.6) 



where is the parity reversal operator. 



The four-vectors l{w) do actually carry a vectorial representation of SL(2, C): 

g>i[w) = i{gw), (2.7) 

where the g in the l.h.s. is understood to be in the vectorial representation of SL(2,C), i.e. to 
be an element of S0+(l,3), the proper orthochronus Lorentz group. Similarly, the vectors w 
carry a vectorial representation of SU{2): 

h>w = hu>, (2.8) 

where the h in the l.h.s. is understood to be in the vectorial representation of SU[2), i.e. to be 
an element of S0(3). Remark that this is the SU[2) subgroup of SL[2,C) which stabilizes the 
time normal f ^ := ( 1 , 0) ^ . 



Let m e be a normalized spinor and m the space component of t(m). Define the 2x2 
complex matrix D[m) := [m,Jm) having m and Jm as columns.^^ It is an element of SU[2). 
Using the previous notation, it can be written as 

SU[2) 3 D[m) := [mjm] = e^i^'-Q-'^'^y (2.9) 



The Livine-Speziale (LS) coherent states are defined as [35] 

\m)j ■-D^i\m)\j,j) = D{m)>\i,j), (2.10) 

where is the usual maximal magnetic number state in the standard j-th SU{2) irreducible 
representation and D'^Xm) is the representation of the matrixD(m) acting on ^^^^.^^These 

11 In vectorial representation D(rn) becomes an S0(3) rotation matrix, written in terms of the Euler angles 
[2^, 29, Zip]. Then, m is just the image of the unit vector z viaD(m): 

rh — D{m) t> z. 

This construction is natural and simple. Indeed, in spin | representation, the spinor associated to the z-axis is just 
(1,0)^ (see footnote 13); therefore for it to be sent into m by a matrix D, the first column of D must be m itself. 
Requiring D e SU{2), gives D = D(m). 
1^ Explicitly: 

{^i^hi=x,y,z being the three generators of SU{2] in representation j. 
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states are coherent in the sense that their expectation value of the angular momentum oper- 
ator ^ is peaked around the vector 7 m e with a relative uncertainty which goes to zero as 
~ 1/^7 the semi-classical large-spin limit. In formulas: 

{m\J\m)j=irh and (m|^^|m)j - |(m|_/'|m)jp =7, (2.11) 

where {m\m')j is the Hermitian scalar product in The LS states can therefore be used to 
described semi-classical three-vectors of length j and direction m e S2. 

Observe that in the previous notation: 

1 1 

m = \m) = \m)i=D{m)>\-,-), (2.12) 
2 2 2 

which holds thanks to\\,\) = [l,QY P 

LS states are an (over-) complete basis of the vector space ^^^^ supporting the representa- 
tion i of SU[2). They allow to write the following resolution of the identity 

^_(2^^+l) dD{m)\m)j{m\j, (2.13) 

JSU{2) 

where the measure appearing in the integral is the normalized Haar measure on SU[2). 

The phase of \m}j plays no role both in the resolution of the identity and in the formula 
associating to \m)j a vector jm. Indeed, its phase has no physical meaning and can be fixed 
to any value. Since the phase of the spinor m translates into a phase 2j(f) for \ m}j, for what 
concerns the definition of the LS states it can be fixed to zero. 
Therefore, we modify the previous notation into: 

\m)j~D[m)>\i,i) with meS2 (2.14) 

and why the resolution of the identity shall rather be written as 

Ij-=(2;-|-l) Am\m)j{m\j-.\ A.m\m) j{m\j, (2.15) 

where dm is the normalized measure on the two-sphere S2, while d^ m := (2j -I- l)dm. When 
no risk of confusion arises, d^- m is simply noted dm. 



Indeed, ||, \) is defined hy \) = \\\, \). Being /i'^"' = one finds ||, |) =(1,0)^. 
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3 BF Theory 

In this section the main steps of this article are briefly reproduced in the simplified context 
of BF theory. This is done, despite their high redundancy, with a twofold purpose: this section 
aims to clarify the ideas lying behind the rather intricate algebra of the calculations in the 
EPRL-FK model, and to justify (or at least strengthen) the reliability of the approximation and 
simplification schemes used. Indeed, such schemes are shown to lead to the correct result 
within BF theory. 



3. 1 BF Graph Amplitudes in the Livine-Speziale representation 



The St/(2)-BF amplitude of a closed 2-complex S^, given an (arbitrary) orientation of its 
edges {e} and faces {/}, is formally defined by 



BF ._ 



e JSU{2) 



where h 



eedf 



(3.1) 



Here Oee5/ stands for the ordered (non-commutative) product over the edge index e along 
the boundary of the face /. Also, egf are the entries of the edge-face incidence matrix: it is 
non-zero only if e e 5/, and it is 1 if the face and edge orientations agree and —1 if they 
do not. Finally, 5 is the delta distribution on the group SU{2] for its Haar measure dh, i.e. 



j^^^^^dh5{hh-')fih) = f{ho). 



The previous expression can be readily generalized in order to include open complexes 
[1]. In this case the boundary of the two complex intersects the external (half-) edges at points 
(nodes {n}) and the open faces at segments (links {I}). The link associated to a certain open 
face / connects the nodes associated to the (half-) edges bounding /. At each link is associated 
an SU[2) group element hi. Therefore, one defines the open complex amplitude by extending 
the previous definition of to open faces. This is done just by formally treating, in the formula 
for h f, the nodes as vertices and the link as edges. Despite of the similarity the interpretation of 
the complex amplitude is different: the hi constitute the boundary data and are not integrated 



over: 



W^^ihi):-- 



n 



dhe 



e JSU{2) J / 



(3.2) 



where the label e applies only to the internal edges, while the label / to both closed and open 
faces. 



Despite of its appealing simplicity, this can be only a formal definition of , since for a 
general complex "ia it contains a priori multiple products of 5(1). One way to regularize these 
divergences is through the Fourier-transform of the St/(2)-BF amplitude 3.2 via the Peter- Weyl 
theorem. In this way it is possible to introduce a sharp cut off at spin j = A. Indeed, the SU{2) 
delta function can be decomposed onto the SU[2) representations labelled by j as 

5ih]=^i2j + l)Trj[h) V/ieS[/(2), (3.3) 

where Trj{h) stands for the trace of the matrix h in representation j, i.e. the 7 -th character of 
h. By using the fact that each of the addends in the previous formula is finite, one can define 
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the regularized version of the S[/(2)-BF amplitude of a complex:^* 



{jf<A] 



nf Y\i2jf + l)Trj^{hf)\, 

e JSU{2) J / I 



(3.4) 



where the sum and the product Y\f are over all faces of the graph, either close or open. For 
future convenience we introduce the amplitudes 



UfKA] I L e JSU{2) J / 



(3.5) 



where the primed summation indicates that the sum is taken only over the close faces, in such 
a way that the ji 's are external data for these amplitudes, while the primed product means that 
the weights (2 ji + 1) are omitted for the external faces: 



/ 



/open 



n^27/+l)Tr/,(^/) 



/close 



(3.6) 



In passing, notice also that the cut-off delta function SaW := X!j<A^^ J ^^^^jW when cal- 
culated at the identity diverges as A^, since Trj(I) = dim^^^' = 27 + 1. 

In order to deal with a formalism closer to the one used in the following sections, it is nec- 
essary to further manipulate Eq. 3.5. First, every group element associated to one edge is de- 
composed into two group elements associated to each half edge: 

he = h~^^^^ht[e)e, (3.7) 

the functions and t{e) denoting the source and the target vertices of the oriented edge e, 
respectively. In this way, there is one group element h^e associated to each vertex-edge couple 
[ve). By construction, h^e always "enters" the vertex. 

Then, for every face, the trace in Eq. 3.5 is decomposed by inserting between each hsf^e)e 
and ht{e]e the resolution of the identity in onto the LS states \ m}j 



Ij = I djm \m)j{m\j. 



(3.8) 



Furthermore, the information about each boundary group element hi is translated into the 
(redundant) information about two LS states m„i and m„'/, being n = s[l) and n' = t[l): 



hi = \mn'i){mni\ + \Jmn'i){Jmni\. 
Hence, the S[/(2)-BF amplitude can finally be written in the LS representation: 



(3.9) 



(3.10) 



i^In the context of BF theory, it is often more effective to use a different, possibly smooth, regularization (e.g. the 
heat kernel regularization, see among others [36, 21]). Nevertheless, here we shall focus on this one, because it is more 
easily implemented in the EPRL-FK model. 
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(3.11) 

Y[ (27/ + 1) Y[ {^e'f\h~]>hye\mef)j^ , 
f "edf 

where the partial amplitudes w^/ have been introduced. Here the labels [ve) and [ve') denote 
the two half edges belonging to the face / ending in the vertex v and such that e' is the edge 
following e according to the orientation of /. Furthermore {m\h\n)j is a shorthand notation 
for the Hermitian scalar product in Remark that in this representation the boundary 

data are more conveniently chosen to be the spins associated to the open faces ji and the 
unit vectors m„i associated to each couple [n, I) composed of an external half edge n and 
an open face I, to which n belongs [n e dl). These boundary data enter the previous ex- 
pression thanks to the following two conventions: (1) if / is the open face such that / e 5/, 
then jf = ji; (2) if / is the open face such that I e df and [ve) (or {ve')) is an external 
half edge ending on the node n ^ dl, then {me'f\h~^^,hve\mef)jf = {me'f\h~l,hue\mni}ji (or 
{rne'f\h~l,hve\rnef)jf = {rri„i\h~l,hve\rri„i)j,,Tespectively). 

In the following, the spins associated to the closed faces of a graph shall be often referred to 
as internal spins. This is because they are associated to faces living in the interior of the graph. 

Finally, remark that in Eq. 3.10 the elements {hj,g} appear always multiplied two by two, in 
such a way that one h,^e multiplies only other group elements at the same vertex. Thanks to 
Haar measure's translation invariance, this fact implies that at each vertex v one of the hyg's, 
say hpeo can always be reabsorbed into the integrations over others h„e- Therefore the inte- 
grations over the various h^eo ^re redundant. This has no consequence for a compact group 
like SU[2). However, it is judicious to gauge fix one group element per vertex in order to 
avoid meaningless divergences in the amplitudes when generalizations to non-compact gauge 
group are envisaged. This will be our case, since the gauge group of the Lorentzian EPRL-FK 
model is the Lorentz group double covering group SL{2, C), which is non-compact. 

3.2 The BF Melon Graph 

The aim of this subsection is to sketch the calculation of the dominant contribution to the 
SU[2)-BF melon graph in the A ^ oo regime, while keeping the boundary data fixed. To do 
this, we shall focus on the limit in which all the internal spins jf are scaled to infinity via a 
parameter A. The idea is to understand the asymptotic behaviour at large spins of the partial 
amplitudes w^g, i.e. of the graph amplitudes at fixed spins. This is the relevant regime for 
studying the convergence properties of the full amplitude: the full amplitude is obtained via a 
sum over six variables, therefore it converges only if the partial amplitudes w^^^ scale to zero 
faster than A~^, and it diverges otherwise. In particular the divergence is logarithmic in the 
cut-offif w^P~A-6. 

It is however important to remark that this procedure is not justifiable on a strict mathe- 
matical basis, since divergences could be hidden in those sectors where one of the internal 
spins is kept small. However, there are at least two independent reasons why to be interested 
in the uniformly large spin regime. First, it is a regime with a sensible geometrical meaning, as 
it shall be shown in the course of this paper. Then, the restriction to this regime is a posteri- 
ori ]ustified by the fact that it is enough to account for the asymptotic behaviour of the melon 



w^Ui'm„i;jf):-- 



nn d 

/ ee3f-'S2 



-Jf 



m 



ef 



nn 



v e:vedeJSU{2) 



dh^e 
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Figure 2: The melon graph. On the right, its faces and extra structure entering the LS represen- 
tation are put into evidence. Triangles point in the direction of the action of the SU{2] elements 
{h, ha,h, ha}- Dots represent insertions of the resolution of the identity. External faces are drawn 
in a dashed line. 



graph in the context of the SC/(2)-BF theory. 
3.2.1 Notation 

The melon graph M (see Fig. 2) has two external half edges, two vertices [v and v), four 
internal edges (labelled by indices a, fo, ••• e {1, . .. ,4}), six internal faces (labelled by unordered 
couples ab), and four external faces (labelled by indices a). Each internal edge belongs to a 
unique external face, whose spin is thus noted ja ■ In the same way the boundary unit vectors 
are noted Ua and fia, depending on the vertex belonging to their boundary. Since the internal 
faces have two edges for boundary, their spins are noted jab = jha- The unit vectors used to 
resolve the identity within the internal edge a, are noted niab if they are inserted in the trace 
relative to the face ab (remark that rhah 7^ rhha), and just rha if they are inserted in the trace 
relative to the external face a. Finally the group elements associated to the internal half edges 
hye and hpe are noted ha or ha, following the obvious identification. The last two group ele- 
ments associated to the external half edges are noted simply h and h. 

With this notation the melon graph amplitude reads 

W^'\ia>na,na)= ^ W^JU a , U a , n a, j ab) (3.12) 

{iab<A} 
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W^^{ia,na,na;iab)~ 



dhdh Yl I dhadha 

_JSU{2f'^ a JSU{2]' 



ab:a<b ^ S 



(3.13) 
drhabdrhb 



5{h)5(h) 



ni n 

Y[ina\h'^ha\ma)jjma\h~^h\na)j 

Yl {Vab + l){mab\h~^hb\mba)j„t,{mba\h^^ha\mab)j„ 



.ab:a<b 



where the 5 functions were introduced in order to implement the gauge fixing conditions dis- 
cussed at the end of the previous subsection (Sect. 3.1). Integrating out h and h, and introduc- 
ing the short-hand notation 



(3.14) 





n dhadha 




Y[ drhabdrhba 




a JsUi2f^ 




.ab:a<bJS2^^ 



and performing the integrals over the rha 's, one gets immediately: 



U^^JUa,na,ha;jab)= \ ^h^tfl 



W{na\hah^^\na)i„ 



Y\ {2iab + l){mab\ha^hh\rnba)jab{^ba\hi,^ha\mab)u 

(3.15) 



ab,a<b 



3.2.2 The BF Large Spin Limit 



In order to give to the right hand side of Eq. 3.15 a form adapt to the large spin analysis, 
we need to write it in such a way that an exponentiated action proportional to the face-spin 
appears. Such a form will very closely remind that of a path integral. This can be easily done 
via the following observation. Since \ j,j) = (^2j l|> |). it is immediate to show that {m'\m)j = 
{m'\m}^} = exp[2j{m'\m}), where \m) = \m)i = m E (see Sect. 2). Therefore, Eq. 3.15 can 
be given the sought form: 



{jab<A} 



W^yAja,na,na;jab) 



w^g{ia,na,ha;iab)~ \ QshOsth 



W{na\hahJ-\na)j, 



W (2;„f, + l)e^«'i 



ab,a<b 



where the face action S^^ was defined 



(3.16) 



(3.17) 



(3.18) 



^^ab ■=2jabMmba\h^^ha\mab)+2jabMmab\h^^hb\mba) 

which is proportional to face-spin jab- 

We also define a total (internal) action S^^ by adding the various internal face actions: 

S«P:= Sfb- (3.19) 



ab, a<b 
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Remark that due to Cauchy-Schwarz's inequality, thie real part of each face action is null or 
negative, and hence also that of the total action 

m[S^^) < 0. (3.20) 

Before delving into the analysis of the semi-classical (large spin) limit of the melon-graph 
partial amplitudes, we shall analyse the symmetries of S^^. First of all, it is clear that at each 
vertex one can left-rotate all the /Za (resp. h^) by an arbitrary S[/(2) elements (resp. ^), without 
affecting Eq. 3.18: 

ha^kha and ha^kha. (3.21) 

Furthermore, one can as well rotate the spinors on a given edge {mab}b,b^a by a fc^ e SU[2), 
provided one simultaneously right-rotates the sets [ha] and {ha} by k~^: 

rriab^kamab ^h,bi^a 

{ha,ha)^{hak-^,hak-^) 

The identification of these symmetries will be fundamental in correctly evaluating the de- 
gree of divergence of the graph. For reasons explained in more detail in Sect. 5.1, these sym- 
metries shall be often referred to as gauge symmetries. 

In order to study the large spin limit of the partial amplitudes w^g, in which all the inter- 
nal spins uniformly scale to infinity, we formally rescale them by a common factor A. In the 
limit A — > 00, it is clear that the integrals are dominated by the stationary points of the total 
action S^p.i^ Following the same steps of [33], adapted at the present case, it is easy to find the 
following stationary point equations: 

Dl(S^^) = iff [ l^^^a = e-^^^^h mab ^3^33 
[ hbmba = e '^"''haruab 

5h,y^ = = 5f^J^^ iff eabjabmab = 0, (3.24) 

b,bj^a 

where ipab and ipab are arbitrary phases, and €ab = -^ba = -f^ab £ {±1} accounts for the 
reciprocal orientation of edges and faces. Finally, the variation of the action with respect to 
the unit vectors thab does not lead to any further equation. 

3.2.3 The Geometrical Interpretation 

The stationary point equations (Eqs. 3.23 and 3.24) can be shown to define (see ??) at each 
vertex a three dimensional tetrahedron in R^. 



'^^Rigorously speaking, these are not stationary points, but rather stationary hypersurfaces on which S^^ is constant. 
However, the contributions of the external faces to the amplitude is not constant on such hypersurfaces (not in every 
direction, at least). Therefore, it will turn out that the partial amplitude will be given at the dominant order by the 
average of the external face contributions on these hypersurfaces. This "averaging" procedure would not have been 
necessary if another gauge fixing were used instead of h — h — I. More on this topic at the end of Sect. 5.1. 

i^More precisely Cat, £ab are a rewriting (adapted to the present notation) of the edge-face incidence matrix (Cg/), 
up to an overall undetermined sign p. See Sect. 8. 
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Each (internal) edge a corresponds to a face of the tetrahedron. Hence one can associate 
to the vertex opposite to this face the same label a. Then the edge lab going from the ver- 
tex a to the vertex b is given by a solution^'' of the stationary point equations via the formula 
^ab '■= ha > eabjabffiab- Remark that by Eq. 3.24 the edges relative to the same face "close" 
forming a triangle, as they should in a geometrical tetrahedron. Moreover, the dihedral angle 
0ab between faces a and b ("around" the edge iab) is, essentially, given by Zeab^at- 

Furthermore, the combinatorics of the melon graph, tells haw the faces of the two tetra- 
hedra have to be identified to form a (closed) three-manifold. This is done along the graph 
edges. In other term the face a of a tetrahedron is identified with the face a of the other one. 
Topologically, the so obtained manifold is a three-sphere.^^ However, subtleties can arise from 
orientation issues. For the time being, it is enough to observe that, since the two tetrahedra 
have the same edge-lengths, they must be congruent or one the parity-reversed of the other. 
Mathematically this choice is reflected at each vertex in the choice of an overall sign (v^) in the 
relation between the angles 0ab and the phases 2eab^ab- 

3.2.4 The BF Amplitude" 

As discussed in the previous section, the stationary point equations 3.23 and 3.24 define 
two geometrical tetrahedra (one per each graph vertex) of edge-lengths given by the internal 

'^''Not any set of spins {jab} admits a solution to the stationary point equations. Explicit conditions for the existence 
of a solution shall be given further. Also, the role of the degenerate sectors of such geometries is not yet completely 
clear. However, see Appendix C. 

i^In order to visualize this fact, it is useful to go one dimension down: glueing two copies of a triangle by their 
congruent sides, and "inflating" the so obtained "sandwich", one clearly obtains a two-sphere. 

i^This section makes use of some results demonstrated only later, in the context of the EPRL-FK model. Most of 
these results apply unaltered to the BF case, just by restricting them to elements of SU{2). However, this shall not 
always be the case. Whenever this happens it will be explicitly remarked in the text. 
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spins {jab}- These tetrahedra can be either the same or be related by a parity transforma- 
tion. In particular it will be shown that the spins [jab] basically determine the niab up to the 
transformation in Eq. 3.22 and, once the niab's are fixed, they also determine up to a parity 
transformation the group elements 

Hab:=h-^hb^SU{2). (3.25) 

This is everything^o is needed to evaluate the melon graph amplitude in the stationary phase 
approximation. 

As a result, the action at the stationary point is (essentially) given by 

Slit = ^^^v + Vu) iabOab, (3.26) 

ab,a<b 

where the 0ab £ [0, ti] are the external dihedral angles associated to the geometrical tetrahe- 
dron of sides {jab}^ and Vi, and Vy are sign factors (±1) assigning an explicit orientation to the 
tetrahedra associated to each vertex of the melon graph. Hence, it is natural to introduce the 
Regge actional for one tetrahedron 

SReggeQ'] := ^ jabOab- (3.27) 

ab,a<b 



Equal Tetrahedra In this case Vy = Vy = ±.\, and the action at the stationary points is given 

by 



nBF 
Stat, eq 



= ±2iS 



Regge 



[J], 



(3.28) 



Moreover, at these stationary points, /z° # h° in any gauge (even if they are simply related), 
where the emphasizes the fact that this group elements are functions of the spins. This can 
be easily shown via the observation that /j^ = ffi^ implies Hab = Hab and therefore (Eq. 3.18) 
S^^ L » = 0, which is not consistent with Eq. 3.28. 

aolstat ' T 



With these elements at hand, it is possible to find a formula for tsymmetrieshe asymptotic 
(large spin) behaviour of the partial amplitudes at the equal-tetrahedra stationary point. 
Consider Eq. 3.17, that we report here: 



W^jAja,na,na;jab)= I ^h^fh 



Y\{na\hah^^\na}j, 



Yl {.2jab + l)e' 



cBF 



ab,a<b 



(3.29) 



In the limit where all the internal spin scale like A ^ oo, this integral (which is on a compact 
domain SU{2Y^ x [SzY^^) can be approximated by stationary phase techniques with 



w 



BF 

,eq± 



n (2;afe + i)' 



ab,a<b 



X-\'^^^'''VF{ja,na,na;iabl 



(3.30) 



2°This is almost true. More complete formula shall be given in the next sections. Our intention now is to skip over 
unimportant technicalities. 

21 This definition is justified by the following fact. Given a closed triangulated 3-manifold M, define for each tetra- 
hedron T of edges {/;} the action SJgggg[/i] := X;=i 'i^Jt'l' hence the total action — ^T'^Regge- the other 
hand, the Regge action for M is defined by S^ggj. — Zj5;, where 5, = 0]^ — 2n. Therefore, S*' = S^gge T^od 2n 
provided that the edges have integer (positive) length. In our case this is almostthe case, since [jab] c |n. 
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where the explicit form of the dependence on the external data is left for the moment unspec- 
ified in the function F, but shall be clarified very soon. Meanwhile, the rest of the formula 
deserves an explanation. 

The first term is a product over the internal faces {ab) of the face weight times the measure 
factors entering the integrals over S2, contributing with one and two copies per face of the 
{"^jab + 1)> respectively. Indeed, recall that dj^^fhat = i^jab + ^)'if^ab> £ind that friab 7^ rni,a- 

In 

A, this factor scales like A^^. 

The second term is the core of the stationary point approximation: for each direction around 
the stationary point in the integration space modulo gauge symmetries, the integral contributes 
by a factor^^ ^ 27r/A, provided that the Hessian of the action 5^S^^ is not degenerate along that 



Now, in order to calculate F in the previous formula, first suppose that the integrand in the 
original expression for w^g were invariant under the action gauge transformations of Eqs. 3.21 
and 3.22. In this case the conclusion would just be that F is the evaluation on the stationary 
phase gauge orbit of the integrand we just mentioned. However, it is immediate to verify that 
the actual integrand is not invariant under the gauge transformation at the vertices (Eq. 3.21), 
because of the external face contribution. The way out of this impasse is to consider the av- 
erage of the integrand on the gauge orbit. At a practical level, this procedure "undoes" the 
previous gauge fixing. Indeed, it result in 



Actually, remark that the integrand would have been invariant, hadn't it been gauge fixed with 
5[h)5[h). Moreover, in the context of SU[2)-B¥ theory, such gauge fixing was superfluous and, 
if it hadn't been performed, Eq. 3.31 would have been obtained straightforwardly.""* 



integration space in Eq. 3.17, and then subtract to it the dimension of the gauge orbit related 
to each solution of the stationary point equations. This must be done, because the action is 
constant on each gauge orbit; hence, the directions spanned by the gauge transformations do 
not contribute to the Gaussian integrals. Therefore: 

rankid^S^^) = 8 x dim(Sf/(2)) -I- 12 x dim(S2) - 4 x dim(S[/(2)) - 2 x dim(SLr(2)) = 30, (3.32) 

eq 

where 8 is the number of ha and ha, and 12 that of rhab over which we are integrating. Mean- 
while, the 4 copies of SU[2) we are subtracting correspond to the freedom of redefining the 
fhab at each edge following Eq. 3.22, and the two more copies are due to the freedom 3.21. 
Remark also that this truly is an upper bound to the rank of the Hessian, which we supposed 
to attain its maximal value compatible with the symmetries. Therefore what we find is a lower 
bound on the degree of divergence. It shall be the comparison of our final result with the well- 
known solution to the problem, to assure a posteriori that such lower bound is attained in the 

^^The 2n's are omitted in our order-of-magnitude treatment. 

2^This conclusion could have been reached by another way, i.e. by performing in the first place a different gauge 
fixing (e.g. via 5{hi)5{h4y), which would break the vertex symmetry of Eq. 3.21, and obtain the integrals over SU{2) 
in the previous formulas from the integrals (no more gauge-fixed) over the group elements h and h associated to the 
external half-edges. 



direction. 




(3.31) 



Before giving the explicit formula for 

eq±' first estimate the rank of the action Hes- 
sian. As already said, to estimate rank[5^S^^), one has to first calculate the dimension of the 
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S[/(2)-BF case. 

Finally, one gets for the partial amplitude: 

"^Xq~^'e±''^''Bge[;l dkdkYl{na\khllh°J-'~k-'\na)u. (3.33) 

JsU{2f^ a 

Unluckily, evaluating the total amplitude Vl/^/eq = S|;„b<A} // eq extremely involved, since 
both the global scale and the dependence on the boundary data is depends on the internal 
spins {jab] in a highly non-trivial way. However, because of the highly oscillatory character of 
exp±2iSRegge[7] for large spins, one can argue that this term gets suppressed. We shall come 
back on this point in Sect. 9.2. 

Parity Related Tetrahedra In this case v„ = -Vy = ±1 and the action evaluated at the 
stationary point is independent on the internal spins and is essentially zero: 

Sstat,p = 0- (3.34) 
Moreover, it is also possible to show, that in this case 

Hab = Hab ^ K = kK' (3.35) 

k being an arbitrary SU[2) element. Therefore, for these stationary points the partial ampli- 
tudes read 

"^i-^M (ikYl{na\ma)ja' (3.36) 

JSU{2) a 

where the fact was used that the estimation of the rank of the action is the same as for two 
equal tetrahedra. This is a very simple expression, and the total amplitude can be easily esti- 
mated. Indeed, the dependence over the internal spins is confined to the global scale of the 
partial amplitude, while has completely dropped out of the terms involving the boundary data. 
Hence: 

K:^= S "^Ip~^'[ dkYl{na\k\na)u, (3.37) 

{jab<A} JSU{2) a 

as a consequence of the summation over the six internal faces. 

At this point, we can put together the two contributions from the equal and parity-reversed 
tetrahedra solutions to the stationary point equations. Accepting the argument according to 
which the imaginary exponential of the Regge action appearing in the equal tetrahedra sector 
causes its suppression, one is lead to write at the dominant order: 

dkYl{na\k\na)j„. (3.38) 

JSU{2) a 

This is the result of our analysis of the SU[2) BP self-energy graph. Notice it is polynomially 
divergent in the cut-off. 



* * * 
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This result can be verified through a simple calculation which is available in the BF case. Start- 
ing from the formal expression in terms of SU[2) delta functions of the S[/(2)-BF amplitude 
(Eq. 3.2), it is indeed possible to first perform all the non-trivial integrals, and then to count 
how many redundant delta functions are left. Only at this point the delta function regulariza- 
tion is introduced (i.e. 5(1) — > 5a(I) ~ A^) and the divergence degree readily calculated. 

In formulas,it is found^^ 



5(I)]M dkYl{na\k\na)j,^A'' I dfc ]^ (3.39) 

JSU{2) a Jsui2) a 

Therefore, Eq. 3.38 correctly captures the dominant behaviour of W^^'"^. 

This is an a posteriori check of our procedure in general and of some of its details in par- 
ticular. E.g. the fact that the rank of the Hessian attains the maximal value allowed by the 
symmetries of the action; or the fact that the equal-tetrahedra stationary point is actually sup- 
pressed; or finally the fact that it is the uniform large spin limit the relevant regime to study 
the most diverging contributions to the melon-graph amplitude.^^ We shall use all these facts 
as a guide in the more complex gravitational case, we shall introduce in the following section 
about the EPRL-FK Spin Foam model. 

To conclude, we add one last comment about S[/(2)-BF theory. It is clear that the method- 
ology employed above is by far not the most economic nor the most effective in this context. 
Despite this, it has the merit of being directly generalizable to the EPRL-FK model, while most 
of the other techniques one can find in the literature make substantial use of the peculiar prop- 
erties of the particularly favourable BF-theory structure. The interested reader can find a much 
more thorough analysis of the melon graph in SU[2)-BF theory, and in particular of the sub- 
leading terms in its amplitude, in [36] . 



Explicitly: 

Y\ 5iH,Hi^')Y\ina\Ha\na). 

ab,a<b a 

It is clear that the delta functions impose Ha — Hb ^ k Va,b. However, this is already implied by the delta functions 
with a — 1 and b free to vary from 2 to 4. The remaining delta functions contribute each with 5(1). There are three of 
these terms, i.e. [ab] e {(23), (24), (34)}. 

2^It could also be mentioned the fact that this formula shows that the degenerate sector, which we did not men- 
tion in the context of BF theory to reduce its dissertation to the minimum, does not contribute at this order. See 
Appendix C. 



n 



a JSU(2) 
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4 Lorentzian EPRL-FK on a 2-Complex 

Starting from this section, we shall treat the physically interesting case, i.e. the Lorentzian 
EPRL-FK amplitudes. In particular this section is dedicated to the general definition of the 
model. 

The cut-off Lorentzian EPRL-FK amplitude Wg of an open two-complex "iS, in the LS rep- 
resentation reads (e.g. [1]) 







nn [ ^ir^^f 




n n [ 


{if<A} 








. V e:vedeJSL{2,C] 



nW) n (^e'f\Yr^g;l,gueYr\mef)jf \ . (4.1) 

where the same notations are used as in the SU[2)-B¥ case (Eq. 3.10). The new ingredients with 
respect to the BE amplitude are the use of SL[2, C) as the gauge group, and the insertion of the 
EPRL-FK -map. This map provides an embedding of the representations Jf^^^ of SU{2] into 
the unitary representation ^(^^-p) = ^iriJ) of SL(2, C): 

Y^ : ^U) ^irJ.j) 
\j,n) ^ \{rj,jli>n) 

where n e {— j, . . . ,7 } is the spin magnetic number. The parameter 7- e M appearing in the def- 
inition of the Y^-map, is the Barbero-Immirzi parameter. The role of the Y^-map is crucial, 
since it is meant to implement the simplicity constraints which, in a Plebanski formulation of 
gravity, transform an SL[2, C)-BF theory into general relativity [3] . 

Remark that in Eq. 4. 1 the face weight is left unspecified. This is because there is no general 
consensus thereabout. Mainly two schools of thought exist: = 2j + 1 [37, 1], as in SU[2)- 
BF theory, or fi[j) = (1 -I- [3, 4], as in a SL[2, C)-BF theory constrained onto the image of 
the Y^-map. We shall keep these two options in mind, while generally noting the asymptotic 
scaling of by 

pt(A)~A'' as A ^00. (4.3) 

Nevertheless, the choice = 2j + 1 seems to us more natural, since it makes the face ampli- 
tude invariant under splitting [37] . 
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^ dm,\m,){m„\ 




Figure 4: The EPRL-FK melon graph. With respect to Fig. 2 lines are added close to the triangles 
indicating group averaging: these lines represent the EPRL-FK Y^-map. Integrations over the CF^ 
spinors {Zab,Zab] are not shown in the figure. 



5 The EPRL-FK Melon Graph: Notation 



In the particular case of the melon graph ^ the notation used here, parallel those intro- 
duced in Sect. 3.2.1: 



Wj({ia,na,na;iab) 



at ria> fiat jab) 

{}ab<A] 

dgdg Yl dgadgfl 

_JSL{2,Cf'^ a JSL{2,Cf'^ 



(5.1) 



Y[ I drhabdmba 

ab:a<b •' Sz 



Y\ drhairialg ^ ga\ma)YrJaima\gJh\na)Y,Ja 
a Js-z 



Yi f^^jab) {mab\gjgb\ mba )y, j„fc {mbalgb^gal mab)Y,.jab 
ib:a<b 

(5.2) 

where we used the shorthand notation {m\g\m)Y^^j := {m\Yj'^ gYj\m) j . 

This amplitude can be set in the form of a path integral. The details of the procedure leading 
to the path-integral formula can be found in [33] (and also [31] , where different - but ultimately 
equivalent - conventions for the splitting of the traces onto coherent states are used). Here, we 
just give the final result: 



a JS2 



Y\ dma{na\ga\ma}Y,.ja{ma\gJ\na}Y,.ja 



ab,a<b 



(5.3) 
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where a shorthand notation for the integrals has been introduced: 



I 



n 

a JSL{2,i 



dgfldga 



n 



ab,a<b S2 



dthabdrhbc 



CP' 



2Jab + l-\ 



n 



(5.4) 

Before defining the new piece of this integration measure, we shall give the new form of the 
face action Sab- 

. , {mba\Zba)^{Zab\mab)^ , • ■ , {^ab\Zab) , 
Sab ■■= Jabln +ir]abin _ ,„ . + 



(Zi,a|Zfoa)(Zai,|Zfli,) 



(Zi,a|Zi,a) 



+7afcln — = = = = \-iYJabi^— = • (5.5) 

{Zab\Zab){Zba\Zba) {Zab\Zab) 

Remark it is still proportional to jab ■ 

The Zab 's are spinors defined by: 

Zab ■= glzab and Zba ■= glzab , (5.6) 
where Zab = Zba ^ CP^ (i.e. a spinor defined up to complex rescaling). Remark thatZab T^Zba- 

Finally, J?^^ is a measure on CP^ (i.e. a scaling invariant measure on C^), defined by 

^ab 



where 



{Zab\Zab){Zba\Zba)' 



2 '^'^ 



(5.7) 



(5.8) 



It is understood that equations analogous to Eqs. 5.6, 5.7, and 5.8, hold at the other vertex 
of^ (tilded quantities). 

The need for the introduction of the CP^ variables {Zab>Zab]> which are absent in the 
S[/(2)-BF case, has to be traced to the necessity of explicitly writing the scalar products in 
the ^^ri.j) = Y^^U) representation of SI(2, C). Indeed, the canonical representation of ^(^^-p) 
is in terms of homogeneous function of two complex variable [31, 33] (see also [38] for the 
general approach). 



5. 1 About the Symmetries of the Internal Action 

It is now important to analyse the symmetries of the total "internal" action:^^ 

ab,a<b 

Such symmetries will be found to be an obvious generalization of those found in the SU{2)- 
BF case. First of all, at each vertex one can left-rotate all the ga (resp. ga) by an arbitrary K e 
SL(2,C) (resp. K), provided one simultaneously transforms the Zab (resp. Zab) appropriately: 

[ JZab^ KJZab [ JZab^ KJZab 

'This action is said to be "internal" because it is related to the internal faces of the graph only. 
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Furthermore, one can rotate the spinors on a given edge {mab]b,bjta by a fc^ e SU{2), provided 
one also right-rotates the sets [ga] and { ga } by the inverse of the same ka : 



mab^kairiab yb.by^a 



(5.11) 



Once more, we stress on how important the identification of these symmetries will be for a 
correct evaluation of the degree of divergence of the graph. 



Since these symmetries are symmetries of the total internal action and of the "internal path 
integral" integration measure 



f^dgadga 



Y[ drhabdniba 



ab,a<b 



n 

ab,a<b 



2jab + 1 



71 



ab ab 



(5.12) 



they shall be often referred to as gauge symmetries. However, this term is inaccurate. Indeed, 
the amplitude of the external faces of the melon graph is invariant only under the transforma- 
tions at the edges"^^ (Eq. 5.11), while it is no? invariant under the transformations at the vertices 
(Eq. 5.10). This is essentially because of the gauge fixing performed on the external (half) edges 
(see end of Sect. 3.1). Therefore, as argued in Sect. 3.2.4 just before Eq. 3.31, the evaluation of 
w^Uab) at the stationary "points" of the action is rather an average on the internal "gauge" 
orbits of the solutions to the stationary point equations. 



If this may sound odd, there is an equivalent way of performing this calculation that does 
not need this "averaging". It consists in gauge fixing to the identity a couple of internal holonomies 
(e.g. gi and g^). In this way, there will be no invariance at the vertices for the internal action, 
and consequently no averaging; nevertheless, the integration over the two SL(2,C) group ele- 
ments associated to the external edges would play the same role as the averaging in the previ- 
ous setting. Also, the degree of divergence of the graph would stay obviously unaltered: there 
would be two less integration on SL(2,C) elements, but also two less SL(2,C) symmetries to 
take into account. Eventually, the two settings are equivalent. 

The reason why we did not choose the latter option, is to keep an higher degree of symmetry 
in the stationary point equations. 



^^Provided also 
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6 Stationary Phase Approximation 

The basic assumption of tiiis section, and of the entire article, is that it is meaningful to 
apply to the EPRL-FK melon graph the same techniques introduced in the previous sections in 
the context of the S[/(2)-BF theory. Therefore, we shall suppose that the uniformly large inter- 
nal spin regime is the pertinent one to be studied in order to understand the dominant term of 
the EPRL-FK self energy. However, one more assumption is needed here. The stationary phase 
approximation can indeed be confidently applied only to functions defined on a compact do- 
main. On the contrary, the EPRL-FK amplitude requires multiple integrals on SL[2, C) which is 
non-compact. Therefore special assumption have to be made on the behaviour of the SL[2, C) 
integrand for large group elements, in order to ensure the applicability of the techniques here 
proposed.^** Supposing all the underlying hypothesis are satisfied, we go through the same 
analysis we sketched in the S[/(2)-BF case. 



By formally rescaling all the internal spins [jab] by a common factor A — > oo, we are lead to 
consider the stationary points of the total action S: 



ab,a<b 

Maximizing29 the real part of the total action (which can be easily shown to be always neg- 
ative or null via Cauchy-Schwartz inequalities) leads to 

3l(S) = iff \ _i^Jf " - (6.2) 

for some phases ip^,(p'^^ [0,27r]. Using Eq. 5.6 and the definition of /, this implies: 

\ If In •- ' (6.3) 

[ gbJmba = ^§^^e'^''^gaJmab 

where ipab '■= Vt~'^a~ ~Vba, and similarly for their tilded counterparts. The vanishing of the 
variation of the action with respect to the variables Zab, Zah, ga and ga gives 

f gbmba=^)Me-'^'"'gamab 
5..,S = = 5,„,S iff . I , (6.4) 

[ gbmba = y^^e "f'^'garUab 

5g^S = = 5g^S iff ^ eabjabmab = 0- (6.5) 

Here, as before, rfiab ■= {mab\d\mab) e S2 and e^j, = -eba = Sab ^ {±1}- The particular 
choice of face orientations we implicitly chose in writing Eq. 5.3 gives [up to a global sign) 
eab = 1 for fl < Z?. Finally, the equations issued by the variations of the niab 's do not lead to any 
further condition. 



The stationary point equations Eqs. 6.2, 6.4 and 6.5 are slightly laborious to work out. For 
this reason, and because of the fact that it is straightforward to apply to our case the techniques 

2^We thank M. Han for having pointed out this fact. 

2^In our conventions there is no minus sign in front of the action within the exponential inside the path integral 
(Eq. 5.3). 
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used in [31, 33], we shall not detail the steps which lead to the stationary point equations. The 
main thing one has to take into account is the fact we are rescaling only the internal spins, and 
not all of them as in the previous references. Concretely, the action to be extremized is that 
of Eq. 6.1, which only involves the data along internal faces of the melon graph. In this sense 
the graph external faces drop out of the "closure" equations 6.5, and no "parallel transport" 
equations-^^ analogous to Eqs. 6.3 and 6.4 has to be considered along such faces. At a practi- 
cal level, it is as if we were treating the large j limit of a closed melon graph with four strands 
only. This is graphically represented in Fig. 4 by the use of different styles for the internal- and 
external-face strands. 

Characterizing the solutions to these equations will be the goal of the following sections. In 
particular, equations 6.2, 6.4 and 6.5 shall be explicitly solved for the variables {fhab}, {ga,ga\ 
and {Zab'Zab} as functions of the spins {jab} only. These solutions shall be finally used to 
calculate the phase of the multiple integral of Eq. 5.3 which defines the melon graph EPRL-FK 
amplitude at its stationary points. 

7 Solving the Stationary Point Equations 

The aim of this section is to explicitly solve the stationary point equations relative to the 
inner faces of the EPRL-FK melon graph (Eqs. 6.2, 6.4 and 6.5) in terms of the spins [jab] only. 
Remark that the quantities relative to the two vertices must satisfy the same equations. As it 
shall be made clear, the sought solution has a neat geometrical interpretation in terms of two 
glued tetrahedra. Each tetrahedron corresponds to a vertex of the melon graph. It is character- 
ized by its edge lengths, which are given by the spins [jab] themselves. For the sake of clarity, 
we anticipate here the result at one vertex, leaving the explicit calculations at each vertex to 
the rest of Sect. 7, and the glueing of the solutions at the two vertices to Sect. 8. 

It is useful to start with the following geometrical remark, whose demonstration can be 
found in Appendix A. Any set of six edge lengths {jab} (strictly) satisfying triangular inequal- 
ities for any subset {jab}b,b^a defines up to global rotations and boosts, and up to space in- 
version, a (non-degenerate) tetrahdron in R^-^. The tetrahedron can be either Euclidean or 
Lorentzian. In the first case the edge lengths satisfy also the condition of positivity of the so- 
called Caley-Menger determinant. 

At each vertex of the melon graph, the solutions to the stationary point equations can then 
be classified into four cases, depending on the particular values of the spins {jabV- 

No Solution If the {jab}b,b^a do not satisfy the triangular inequalities at each edge a, then 
there is no solution to the stationary point equations (in particular to the closure equations). 
The corresponding graph amplitude results exponentially suppressed with the spin scale A. 

Euclidean Sector If the triangular inequalities are strictly satisfied at every edge, and the 
Caley-Menger determinant of the spins [jab] is positive, then there are (at each vertex) exactly 
four solutions to the stationary point equations equations, which correspond to two different 
geometries characterised by a sign choice. These two geometries can be interpreted as the two 

^"As it shall be made clear later, "closure" equations are equations like Eq. 6.5, while "parallel transport" equations 
are equations like Eqs. 6.3 and 6.4. 
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different, parity-related^^ Eculidean tetrahedra one can build out of the six side lengths {jab}- 
The two parities are characterised by a sign Vy = ±1 (Eq. 7.21). 

Lorentzian Sector Similarly, if the triangular inequalities are strictly satisfied at every edge, 
but the Caley-Menger determinant of the spins [jab] is negative, then there are (at each ver- 
tex) exactly four solutions to the stationary point equations equations, which still correspond 
to two different geometries characterised by a sign choice. However, these two geometries can 
now be interpreted as the two different, time-reversal-related "Lorentzian tetrahedra" one can 
build out of the six (Minkowskian) side lengths {jab}- Remark that > 0, hence the recon- 
structed tetrahedron has space-like sides and faces. Once again the two possible geometries 
are characterised by a sign Vy = ±1 (Eq. 7.21). 

Degenerate Sector If the triangular inequalities are satisfied at each edge by the {jab} b,b=iia , 
but they are saturated at least at one edge, then the graph is said to belong to the (geometri- 
cally) degenerate sector. It shall as well be called degenerate, the sector spanned by the {jab} 
strictly satisfying the triangular inequalities, whose Caley-Menger determinant is null (zero 
volume tetrahedra). The study of this sector is left for future work, however see Appendix C. 

* * * 

This section is organized as follows. In the first part we shall completely solve the full set of sta- 
tionary point equations in a purely algebraic way. In this way the existence of the two different 
families of solutions is shown (Euclidean, Lorentzian). In the second part, we shall provide a 
geometrical interpretation to the stationary point equations for the non-degenerate Euclidean 
and Lorentzian solutions. Such an interpretation shall allow us to give an explicit meaning to 
the different quantities involved in the solutions to the stationary point equations at the end 
of this section. 

All along the calculation we shall take advantage of the possibility of solving all the equa- 
tions explicitly, trying to add a flavour of concreteness to the more abstract works of [3 1, 33] . 

7.1 Algebraic Solution 

At an edge a, Eq. 6.5, hereafter called the "closure equation", implies that the three vectors 
{fnab}b,b=^a coplauar. Therefore, we can use the freedom of SU[2) rotating the spinors at 
each edge (see Eq. 5.11), in order to gauge fix the plane the vectors lie on, to be the xz-plane. 
In the notation of Sect. 2, this means that the niab's take the form^^ 

m.. = e'^-f V (7.1) 

^^This can be seen as the change in sign of the z-axis, corresponding to 17 —r]. 

32Remark that this is the case for any choice of ipab — ^k, k& {0, 1,2,3}, possibly at the cost of redefining the phase 
(pab and the angle 6ab- Explicitly: 



4>ab^(t>ab + -k, and dab^[-'^f dab- 
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Remark that the mat's are normalized spinors, according to what is needed for the resolution 
of the identity Ij = djm\m)j{m\j (see Eq. 2.15). 

Via a further rotation of the {rhab]b,b=iia within the xz-plane, it is possible to set Bat = 
for a particular b. Then, it is straightforward to solve the closure equation for the differences 
9^^ := [dab — Sac)> provided that the {jab]b,b^a respect the triangular inequalities (otherwise 
no solution to this equation exists). This gives 

cos20,«^ = -eatea /"'^^'^'.'-''"' ; (7.2) 

^Jabjac 

furthermore, the closure equation fixes the relative signs between the previous differences in 
the following way: 

sgn[eafeSin26'^^J] = sgn[e„dsin26'^^^]. (7.3) 

where in the last two equations the convention is used, that different letters stand for different 
values of the indices. This convention should be used throughout the rest of this section. 

It is now convenient to define a new set of angles {ifab}, in such a way to reabsorb the 
dependence from the {eat}'- 

^ab-Oab + ^^n ^ cosZf^^ = e abe ac cosZOi;^, sm2^^^ = eabeacsm29l^^ (7.4) 
This new definition allows to write the previous two equations in the cleaner form: 

cos 2^^^ = - ^^b+Jlc-jad (7 5^ 

^Jabjac 

sgn[sin2#^J=sgn[sin2#«^]. (7.6) 

The first equation has two solutions for 2#^^ (modulo 2n) related by a sign change.^^ The sec- 
ond equation fixes the relative signs of these angles at a given edge {a, in this case). Therefore, 
at fixed {€ab}, exactly two solutions per edge exist. 

The attentive reader has surely noticed that the previous equation simply follows from the 
interpretation of the vectors {eabiabniab\b,b^a in terms of the sides of a triangle, as dictated 
by the closure equation. Nevertheless, we shall postpone the geometrical discussion. 

Now that the closure equations have been solved by assigning a particular value to the 
mab's, we can focus on the other two stationary point equations. These shall allow us to solve 
for the SL(2,C) group elements {ga}- To do this, we rewrite equations Eqs. 6.3 and 6.4 as a 
unique matrix equation; 

(ll^fcnllg-iyab A 

^ WZbar J 

where D[m) := [m,Jm) = e^^'^'^e'^^^ye^'l'^^ is the SU[2) matrix with the two (normalized) 
spinors m and Jm as columns.-^* Remark that in this formula the niab are gauge fixed in such 
a way that j/'ab = 0. It is then convenient to define 

WZabW 

r]ab^-r]ba ■=i^ab-(pab + 4>ba) + in eC, (7.8) 



^^While it has four solutions for 
34Cf footnote 11. 
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and 

Gba ^ ■- g-'ga eSL(2,C) . (7.9) 

In the next section, it shall be shown that these two quantities have a neat geometrical inter- 
pretation. 

In terms of the new variables, Eq. 7.7 simply fixes the form of the Gab'^'- 

However, most of the physical information lies now in the fact that the Gab = ga^gb are not 
all independent from one another. Indeed they must satisfy three consistency equations: 

■ 

GabGbcGca = I 
- GacGcdGda=l ■ (7.11) 
GabGbdGda =1 

V 

A fourth equation follows from the previous three: GbcGcdGdb = I- 

An important feature of these equations is that they involve only the differences between 
the angles at a same edge of the graph. E.g., the first one involves the differences Q^^ = 
{Gba — Qbc)y ^ba ^'^^ ^cb' already observed, such differences are essentially determined by 
the spins jab's via the closure equation. This means that Eqs. 7.11 can be considered as equa- 
tions for the complex variables rjab- 

Also, notice that since the niab's are in common between the equations relative to the two 
vertices, 3[fiab) = 'fab - (pat + 4>ab, for the same (pab and (pba appearing in 3[riab)- This fact 
shall be of importance later, where it will let the unphysical phases 4>ab drop out of the final 
formula. 



Now, in order to solve this system in Eq. 7.11, we start by analysing one single equation, e.g. 
the first one. Using Eq. 7.10, it can be cast into the form: 

gi7fcflO-,gie*^,0-,,g)7,fcC7;: ^ Q-ie;^l,C^y Q-nacC^Z Q-ie;;^C7y ^ ^y_j2) 

where the differences between the angles 6^^ = [Qba — dbc) are put into evidence. 

One way to solve this equation is to explicitly write the two matrices appearing on its left 
and right hand side 




e I"': cos9",cos9f —e'^"'^^sm6",sin0!: — ei^^ sin0'icos05 — e cos01sin0,'^ A 

CO ba cb ba cb ba cb ba 1 

e'l"' cos 9^^ sin 6^^ + e"''"' sin 9^,^ cos 9^^ e''«' cos 9^,^ cos 9^^ - e""'"' sin 9^^ sin 9^^ J ' 



and taking the sums and the differences of the elements along the two diagonals: 

(cH-r]ba - r]cb)cos{-9^^) = chr]ac cos{9^^ + 9^^) 
ch{-r]ba + r]cb)sm{-9^^) = chrjac sinl^^"^ + 
sH-riba - r]cb)cos{-9^J = shr]ac cos(0^'^^ - 0^"^) 
sh{-r]ba + r]cb)sin{-9^J = shr]ac sin{9l^^ - 9^^) 



(7.13) 
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This system of equations can now be easily solved for chlrjac, yielding 

cos20.^ cos20^. -cos20^ 

ch2riac = ".^^^ ^. (7.15) 

sm20^^^sin20;^ 

And using Eq. 7.4: 

cos2#^ cos2#'?, +cos2#^^ 
sin2#^^sin2i^«j, 

Here we defined the function Aadj] e K, which depends only on (all) the spins {jab} (as well 
as on a choice of sign at the two edges a and c. See the discussion after Eq. 7.4). In particular 
it does not depend on the choice of the {eab}- 

Since the function A[j] is real, 2r] can take only particular values (modulo 2i7i): 

[Eucl] ifz\e [-1;1], then 217 = ±i acosZ\, 

[Lorl] ifz\E [l;oo],then2)7 = ±achzi, 

[Lor2] if A e [-00; -I], then 2/? = ±(achz\ + i7r), 

where we neglected cases in which particular numerical coincidences among the values of the 
[jab] are such that A e {±l,±oo}.^^ Remark that the number of solutions for 17 (modulo 2i7r) 
itself doubles with respect to the number of solutions for 217. 

Solutions of type [Eucl] ([Lor]) are said "Euclidean" ("Lorentzian") and will be shown to cor- 
respond to the Euclidean (Lorentzian) geometrical sector. Indeed, the parameters 217 will be 
shown to be related to the dihedral angles among faces of a tetrahedron: real and imaginary 
values of this angles correspond to pure rotations and boosts, respectively. 

Now that Eq. 7.16 gives the norm of every 2r]ab as a function of the {jab}: one would like 
to control also its sign. First of all, notice that if {Zrjab} solves Eqs. 7.14, then so does {—2rjab}- 
Therefore, there is no way to extract from these equations the absolute sign of the 2r]ab's. Nev- 
ertheless, it is possible to extract their relative signs. 

A preliminary step in this direction is to show that if one of the rjab'^ is in the Euclidean 
sector, then so are aU the others. To do this, it is enough to add among them the products of 
the first by the fourth and the second by the third of Eqs. 7. 14, to obtain 

sh2r]ba _ sh2r]ac 



or, equivalently: 



sin20,^„ sin20,^, 



sh2riab sh2r]c 



eab . „cic . „„^ • (7.18) 

sm2T?^j, sm2#^^^ 

As an immediate consequence of this equation, it is found that if 2rjab £ iK at one face {ab) 
(case [Eucl] in the list above), then so will be at any other face. In this case Gab ^ Sf/(2) Vflfo 
(Eq. 7.10). Otherwise, the Gab'^ are genuine elements of SL(2,C) (even if not the most general, 



^^These values correspond to the degenerate geometrical sector. As it shall become clear later, A — ±.\ implies a null 
volume of the tetrahedron with no degenerate faces, whUe A — ±00 implies a null volume of the tetrahedron because 
of the degeneracy of some of its faces. 
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Eq. 7.10 and cases [Lorl] and [Lor2] of the list above). 

To proceed with the analysis of the signs, it is first necessary to perform another gauge 
fixing on the {niab}- The will is to fix the signs of the sin2i?^^'s, which appear in Eqs. 7.6, 7.16 
and 7.18, in some globalway. This shall be geometrically interpreted as the fixing of a coherent 
orientation throughout the tetrahedron. 

We start by noticing that at one edge a, the signof all the (see Eq. 7.6) can be easily changed 
by acting on the {niabh.b^a with ka = e't°'^ e SU{2y?^ 

m.. = e'^- [ ] e^^--'f [ '^'^'^-'f'] ] ^h + a. (7.19) 

\^ ^va.Qab J y sm{-9ab) J 

Indeed, this implies (modulo 2ti) : 

2^ab --^ -'itab ^b^a and therefore 2#^^ -2#^^ \lb, c^a. (7.20) 

The gauge fixing then goes as follows. Using the completely antisymmetric tensor Cabcd such 
that ei234 = +1, and an assignement of a (different) number from 1 to 4 to each edge, the gauge 
transformation Eq. 7.19 is performed at the edge a if and only if sgn[sin2#J^^] ^ eabcd, where 
d is the only element of the set {1, 2, 3, 4} \{a,b, c}. Remark that this prescription is consistent 
with Eq. 7.6 (which was in turn implied by the closure condition) since eabcd = ^acdb- 

Within this gauge choice, Eq. 7.18 acquires a very simple meaning when combined with the 
explicit solution for chZrjab (Eq. 7.16): 

2eab'nab = Vv0ab{j) modZin, where ©afoCj) e | J^^'^^^^ .^^ ^^^^^j (7.21) 

i.e. where 0abij) ^ i[0, tt] in the Euclidean sector ([Eucl]), while Oabij) ^ + {0>i^l) in the 
Lorentzian one ([Lorl] or [Lor2]). Moreover, Oabij) is a function of the spins only, and in partic- 
ular it is independent of the eab's- Finally Vy e {±1} is an overall sign that can be independently 
chosen at each vertex {v,v]. 

The last step in the solution of the stationary point equations, is to find the solutions for 
the rjab's, instead of the Zrjab as in Eq. 7.21. This step is not a completely trivial, since the 
rjab'^ are defined modulo 2i7r. Therefore, the possible solutions for rjab are both | x [2rjab) and 
|x(2?7ai,) + i7r. However, Eqs. 7.11 put some constraints on these choices, which are not all in- 
dependent from one another. In other words, the choice at one face [ab) influences the choice 
at the other faces. Unluckily, there is no geometrical interpretation of this choice. Indeed, as 
it shall be clear from the discussion in the next section, the geometrical interpretation stems 
from the spin 1 representation of the parallel transport equations, which need on the contrary 
be solved in their fundamental (spin |) representation. 

Unfortunately, the way the transformation rjab rjab+i'^ influence the other {ricd}{cd)^{ab] 
is not unique, and various solutions are possible. Indeed, it is not hard to convince oneself via 
Eq. 7.14, that if e.g. rju rjiz + in then one of the following cases is given: 

• also ri23, J734, rju acquire an m; 



^SRemark that this gauge fixing is compatible with the previous one, done at the very beginning of this section. 
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• also ri23, J724 acquire an in; 

• also 57 13, )7i4 acquire an in; 

• also )7i3, (724, J734 acquire an in. 

We shall not go into the details of these solutions, which may be cumbersome to write down 
explicitly. However, we shall generally note the fact that "some" Iti's may be needed "here and 
there" with the following symbolical expression 

r]ab = ^Vueab0abU)+[in]ab mod2i7i, where 0afo(;) e | .^^ ^^^^^j (7.22) 

Luckily, in Sect. 8 it will be shown that the detailed knowledge of these terms is not necessary. 
This is because of some global properties of the graph amplitude, which will allow us to bypass 
this issue. 

7.2 Geometric Interpretation 

Since the work of Ponzano and Regge [2] , the standard idea behind the large spin limit of 
spin foams is that of recovering a classical discrete geometry. We shall show that also in the 
case of the EPRL-FK melon graph with fixed external spins, such an interpretation is possible. 
Since our analysis practically concerns only three out of four strands per edge, it shall be re- 
lated only to a three dimensional geometry, rather than to a four dimensional one as expected 
for more complex graphs [33] . 

At the beginning of the previous section, it was already mentioned how the closure equa- 
tions (Eq. 6.5) define at each edge a geometrical triangle a embedded in M^, and how the latter 
is determined by its side lengths {jab}b,b^a up to rotations. It was then algebraically shown 
that it was possible to gauge fix such triangles to lie in the -plane. In the following, we shall 
make explicit how the other stationary point equation can be also interpreted in terms of these 
triangles: they will be seen to encode their parallel transports from the "edge reference frame" 
to the "vertex reference frame", in such a way to constitute in the latter frame the faces of a geo- 
metrical tetrahedron. Such a tetrahedron shall be Euclidean (embeddable in M^) or Lorentzian 
(embeddable in M^.i)^ according to the values of their edge lengths (or, equivalently, to the val- 
ues of the zi[7]'s in Eq. 7.16). 

On a concrete level, the goal of this section is to give a neat geometrical meaning to the 
rather abstract functions Oabij) of the previous section, which will be found to be anything 
but the dihedral angle (respectively boost parameters, in the Lorentzian sector) among the 
faces of the tetrahedron associated to each graph vertex. During this quest, the following geo- 
metrical meaning will be found for the Gab's- their spin 1 representation will be identified with 
the rotations (respectively Lorentz transformations) taking from one face of the tetrahedron to 
another. 

Despite their algebra and interpretation are very similar, the Euclidean and the Lorentzian 
sectors shall be treated separately since the beginning, because of technical reasons. 
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Figure 5: A graphical representation of the vectors {ha> lab} in the Euclidean sector. Remark how 
they close into triangular faces (Eq. 7.25) and into a tetrahedron (Eq. 7.26). 



7.3 The Euclidean Sector 

In the Euclidean sector, %[riab) = 0, and therefore Gab = Sa^Sb ^ SU{2). Hence, the two 
SL(2,C) parallel transport equations (Eqs. 6.3 and 6.4) become equivalent to one another. 
Moreover, such a solution for the Gab s implies the following form for the ga 's: 




ga = gh[ 



Vfl, with hi E SUiZ), and any g e SL(2, C), 



(7.23) 



where the must be understood to be determined by the jab's. 



Define now 




(7.24) 



Remark that jab = jba, ^ab = —^bay while niab has a priori no relation with rriba other than 
Eqs. 6.3 and 6.4. Also, the symbol eab = —^ba is the same which appears in the closure equa- 
tion (Eq. 6.5). 



The closure equation (Eq. 6.5) being written at each edge, can be acted upon with ha e 



SU[2): 




(7.25) 



b,b^a 



Also, in the Euclidean case, Eqs. 6.3 and 6.4 when written in spin 1 representation, read: 



hl>Iab = -hl>Ib^ 



(7.26) 



In the last two equations, h^ , h° 



are understood to be elements of S0(3). 



Equation 7.25 implies that the three vectors > iab} define, at fixed a, the sides of a tri- 
angle; in turn, Eq. 7.26 implies that they are moreover coherently identified among them, so to 
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Figure 6: The figure represents the action of a SLA(2) element h = exp rj'az on the dashed triangle 
lying in the xz plane. A second (fixed) triangle is also shown. 



form an (oriented) tetrahedron.^^ The latter identification is due to the combinatorics of the 
spin foam vertex (see Fig. 3). 

To extract the geometric content of the parallel transport equations, i.e. in order to cal- 
culate the dihedral angles between the faces of the tetrahedron, define the oriented normal 
vectors to the triangles 

^ 1 ^ ^ 

Fa~--eabcd^ab^tac. (7.27) 

Here no sum on repeated indices is performed, and eabcd is the completely antisymmetric 
symbol used for the gauge fixing of the previous section. Notice, that this gauge fixing, through 
the previous definition, also fixes an orientation for the tetrahedron and its faces, as it was an- 
ticipated in the previous subsection.-^^ Furthermore, notice that this expression is independent 
of the choice of the two sides of the triangle a involved, thanks to the closure equation and to 
the antisymmetry of the eabcd- Now, the cosine of the internal dihedral angle 0^^^ e [0, n] be- 
tween two faces of the tetrahedron, is given by minus the (normalized) scalar product of the 
two outwards (or inwards) pointing normals to the faces of the tetrahedron, which are e.g. 
ga > Pa and -gb>Fb, therefore: 

cose',:="'°"^°'"'°"^'' = ^-'°°'"^''. (7.28) 
where Eq. 7.23 was used. 



^^For example, one possible realization of such a tetrahedron is given by the following positions of its vertices: 
X\ — (the origin of R^), £2 =£34, X3 =£42 and £4 =£23- Here we named the vertices with the index of the face opposite 
to it, and orientations are chosen as in ??. 

^^Indeed, three vectors satisfying the closure equations, can be organized in two different ways, each of which 
provides a different orientation for the triangle. Therefore, asking for Fa to be the oriented normal to the triangle, 
picks out only one of these two possibilities. 
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By gauge fixing the triangles a and b to lie in the xz plane, we have 

Pa=Aay, where Aa:=-^eabcdy-(/abX(ac) = ^eabcdSm2t'^^>0, (7.29) 

thanks to the gauge fixing condition. Also (in spin 1 representation)^^ 

G^t = e-^^"'^y e^<b^' e^^'^'^y (7.30) 

where we used the fact that rjat = 3(i?afc) = '7^^- Hence, a straightforward calculation of 
Eq. 7.28 gives 

^'^^^ab = ^°^^^'ab = C0S2((/?ai, - (pab + 4>ba)- (7.31) 

This formula, combined with Eq. 7.21, makes the connection among the parameters 2r]ab's 
and the tetrahedron internal dihedral angles 0^^, in the Euclidean case: 

0abU) = i0ab [Eucl]. (7.32) 

At this point, the reader may be worried about the appearance of the non-physical phases 
(pab in the previous formula. However, we shall notice that such phases make no appearance 
in any formula with true physical significance (e.g. a transition amplitude), but just in those 
formula which constitute an interpretation of the transition amplitudes (see the second step 
in Eq. 9.2). 



7.4 The Lorentzian Sector 

The structure of this section is exactly the same as that of the previous one. They differ 
just because of some technicalities. In particular, what here is called tab will not represent any 
more the side of one of the triangular faces, but rather the unique normal to it on the plane 
defined by the face itself. Since this fact implies minor changes both in the algebra and in its 
geometrical interpretation, the decision was made to maintain the same symbol for the two of 
them. 



In the Lorentzian sector, Gab = ga^Sb ^ SL[2,C). Equations 6.3 and 6.4 implies that the 
ga 's have the following form: 

= gg° Vfl, with g° e SL{2, C), and any g e SI(2, C). (7.33) 

where the g^, as much as the Gab> must be understood to be determined by the jab'^- 



^Recall the formulas for the Lorentz generators (spin 1 SL(2,C) generators): 
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Now, for every couple {ab), define the four-vector^" 

(■ah~^eabiab*i[My)M[mab)M{Jmab), whcrc My := ^ I j =>t(My)= 1^ ^ 

(7.34) 

where it shall be understood that the niab'^ have been gauge fixed as in Eq. 7.1. Moreover, in 
the previous equation the wedge product stands for the anti-symmetrized tensor product, and 
the asterisk for Hodge duality.*^ In coordinates, this can be written as 

t'ab = ^abjab ^ ^' JKlt' m'ab' (7.35) 
J,K,L 

where the following symbols have been introduced: 

,'=(„i), /=(;) and ™;.=(;^J. 

Therefore, it is clear that f^fo is a spacelike four- vector, and that it lies in the three dimensional 
hyperplane orthogonal to . Here it forms an orthonormal triad together with y and ifiab- 
Also, its length is \\£ab\\^ \f(\J^i = jab- 

Then, we wedge the spin one representation of the stationary point equations in Eqs. 6.3 
and 6.4, with each other and with the four vector t(My). This simply gives, in terms of the lab- 

gl>iab = -gl>eba, (7.37) 

where we used the shorthand notation g > l[z)Ai[z')Al[z") := [g®g®g] > L[z)AL[z')/\i[z") = 

iigZ)AL{gz')AL{gz"). 



Furthermore, from the closure equation (Eq. 6.5) for the fhab can be immediately deduced 
a closure equation for the tab (this can be most easily done via Eq. 7.35). On this equation one 
can then act with the group elements g° , to obtain 

X g°>^«f'=0- (7.38) 

b,bj^a 

The natural geometrical interpretation of Eq. 7.38 is that the {g^ > iab}b,b^a define in E^'^ 
(up to translations) an oriented triangle a with sides g° > lab- In the same terms, Eq. 7.37 
shows that the sides of these triangles are coherently identified in such a way to form an ori- 
ented tetrahedron. 



''"We allowed ourselves an abuse of notation here, since we are assigning here the symbol iab to a different entity 
with respect to the previous section. However, in doing this, we want to stress on the analogous physical interpretation 
of the two quantities. 

''ipor definiteness: 

ViA---AV„ = — Sign((7)l'cr(l) A ■ ■ ■ A Va{^n) , 

and 
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Figure 7: The figure represents the action of an SL{2, C) element g = exp rj^cTz on the dashed 
triangle lying in the xz plane. A second (fixed) triangle is also shown. Remark that with respect to 
Fig. 7 the axis on the plane are rotated. 



As in the previous section, in order to extract tiie geometrical content of the parallel trans- 
port equations, we define the face bivectors as*^ 

Fa~-^eabcd*iabA£ac- (7.39) 

Using Eq. 7.35, this can be recast in the form 

Fa=2AatAy, where Aa :=-^eabcdy-(.iab x£ac)>0. (7.40) 
Or equivalently, in covariant form 

Fa=AalUMy)M{My). (7.41) 

In turn, in spin one representation: 

G^i, = e-^e,bfy Q2v,l,-y^. Q2e,,fy _ (7.42) 

where rj^ and rj^ are the real and imaginary parts of e C, respectively and where we used the 
fact that in the Lorentzian sector 17^^ e {0, n}. 



Thus, the scalar product between two face bivectors is: 

(gg > FaYKgb > Fb)lJ ^ Fl^jGab > Fb)lJ 

WFaWWFbW WFawmw 



chZr^f,. (7.43) 



Analogously to the Euclidean case, this formula gives the dihedral hyperbolic angle between 
two faces of the Lorentzian tetrahedron:*-^ 

ch0^, = ch2r?f,. (7.44) 



''2 As before for the £ab% here we commit an abuse of notation. 

''^Since all the face bivectors are future pointing (Eq. 7.40), we are always in the thick wedge case of [31]. 
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Therefore if 0^^ is the absolute value of the boost parameter among the faces a and b (internal 
hyperbolic dihedral angle), using Eq. 7.21 we get 

0abU) = 0'ab [Lor]. (7.45) 



8 Graph Face Orientation 



Before estimating the radiative corrections to the EPRL-FK propagator, we think a brief 
comment about face orientation is in order at this point. 

The EPRL-FK model is invariant under face orientation reversal of the internal faces (and 
in some sense "covariant" with respect to the face reversal of the external faces) . This fact had 
already been emphasized in [39] , nevertheless we shall give a brief demonstration of this fact 
in Appendix B. 

Furthermore, if the coherent states are eliminated from Eq. 4.1, using the formula for the 
resolution of the identity,^^ then even the integrand of each face amplitude is invariant under 
face reversal (see Appendix B): 



vedf 



is invariant imder / — » / ^ . (8.1) 



However, as soon as the coherent states are introduced, a face orientation is automatically 
picked out, and the integrand is not any more face-reversal invariant. 

Once the real part of the face action 3l(S /) is set to zero by one of the stationary point equa- 
tions, and the CP^ spinors associated to each face-vertex couple (z^y in the notation of [33]) 
can be expressed as functions of the coherent states rrief, it is immediate to see that 

Sf-. [z,f, m,f, g,e] \z^^^^„^^ = - Sf[z.f, m,f, g,,] l^^^^^^^^ . (8.2) 

At the same time, the closure equations themselves are not invariant under face reversal, 
since eef[v) —€ef[v) [iv, e e d f when / /^^ As a consequence the g~^,gve = G^, [v) 
solving the parallel transport equations depend on the face orientation as well. As it shall 
become clear from the explicit example of the melon graph, this two changes compensate each 
other, giving - as expected - a face-reversal invariant phase at the stationary points. Indeed, 
the way we wrote the melon graph amplitude in Eq. 5.2 automatically implies that the faces 
are oriented counter-clockwise and that, once a global choice of v e {±1} has been made: 

e,, = -e,, = | _^ .^^^^ . (8.3) 

9 Radiative Correction 

In this section, we shall use the analysis of the saddle point equations previously per- 
formed, to evaluate the melon graph at its dominant order. In particular, we shall evaluate 

'%^^,^d.m\m)j{m\j. 



38 SELF-ENERGY OF THE LORENTZIAN EPRL-FK SPIN FOAM MODEL OF QUANTUM GRAVITY 



the partial amplitudes at fixed (large) spins in Eq. 5.2 at their stationary points, as well as es- 
timate the rank of the action Hessian at these points so to understand the partial amplitude 
scaling at large spins. Finally, we shall argue how the sum over all the possible internal spins 
affects the final result. 

Formally, the argument goes as follows. First, we rewrite the truncated amplitude as a sum 
of partial amplitude at fixed spins: 

WgUa,na,na)-. ^ a, fla, na', j ah) ■ O-D 

{ia<A} 

Then, we claim that the procedure studied in Sect. 3 in the context of BF theory applies also to 
the EPRL-FK amplitude. Summarizing, since the face weight ^[jab) scales with some positive 
power > 1 of the spin, we claim that the relevant contributions to the sum in Eq. 9.1 come 
from the configurations where the internal spins are all large. The formal way to achieve this, 
is to scale the internal spins jab by a factor A ^ oo, therefore "picking" the stationary points of 
wj/ as the relevant contributions. 

In the next subsections we shall evaluate the partial amplitudes wj{ at their (geometrically 
non-degenerate) stationary points. As it was shown in the previous section, the geometry type 
of the stationary points (Euclidean or Lorentzian) depends on the value of the internal spins 
jab only. Therefore, at the two vertices of the melon graph, the two geometries are of the same 
type. 

9. 1 Action at the Stationary Points 

In order to calculate the action at its stationary points, it is enough to combine together the 
stationary point equations, their solution (Eq. 7.22), as well as the result of the previous section 
about the graph face orientation. 

We first substitute into the total action formula (Eqs. 5.5 and 5.9) the solutions to the sta- 
tionary point equations Eqs. 6.2 and 6.4 (essentially the parallel transport equations), and their 
parametrization in terms of the complex parameters rjah ^ C introduced in Eq. 7.8: 

Sstat = i jab 

ab,a<b 

= i ^ 2jab[-%riab-fjab) + r^{r]ab-'nab)\- (9-2) 

ab,a<b 

Notice that the formula was used %rjab — i)ab) = fab — <^afo, which stems from the fact that the 
same phase (j)ha — (pab appears in both rjab and fjab- therefore the arbitrary phases 4>ab's play 
no role in the relevant formulas. 

At this point, one just need to combine Eq. 7.22 and the result about the face orientations 
of Eq. 8.3 and to insert them in the previous formula, to first obtain 

1 1 

r]ab = -VvV 9ab + [in]ab and flab = --Vvv 0ab + [in]ab vjhena<b, (9.3) 

and then 

Sstat = -VK+Vs)i jah[%0ab) + r^{0ab)]-V ^ Zjabiin]^^,. (9.4) 

ab,a<b ab,a<b 



- f, WZabW , \\Zba\\\ 

-2[fab- fab) + 2r\ In— --In j-^— 

V WZbaW \\Zab\\ J 
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In the last equation we made use of the fact that 2 jab e N to introduced the new symbol^^ 

[i^]afo-=[i^]afe + [i^lafc mod2i7r. (9.5) 

Finally, remark that the choice of the sign —v is irrelevant, since it can be absorbed in that 
of the signs and Vp . Therefore: 

Sstat = i(v, + Vs) jab[%0ab) + r^{0ab)]+ ^ ^J^bMlb' (9-6) 

ab,a<b ab,a<b 

Notice that this expression is manifestly invariant under face-orientation reversal. 



9.2 Euclidean Stationary Points 

In the Euclidean case, irjab^^> and 0ab = ^^ab • There are essentially three possible phases 
for the stationary points, depending on the choice of Vy and Vy : 

1. V, = -Vo = ±1 ^ Sstat = + llab,a<b '^Ub[m]li,-; 

2&3. y, = V, = ±1 =^ Sstat = ±i:ab,a<b Vab {iO^^b + i^Kb) = ±2iSigge[7]+i:„b,„</, 27„,[i7r]«,. 
Here S^^^^^[jab] stands for the (classical) Regge action for an Euclidean tetrahedron of sides 

[jab]- 

For the stationary points of cases 2 and 3, the sum over the spins which is needed to com- 
pute the total amplitude can be viewed as a discrete path integral weighted by the Regge action. 
In the large spin limit where all the jab 's scale as A — > oo, the exponential of the action is rapidly 
oscillating, and the path integral will be dominated by its stationary points, i.e. by the classi- 
cal motions of the Regge action for a tetrahedron. In this case, the well-known equations of 
motions of the Regge action read ©f^Q) = Vafo [40], which have clearly no solution in the 
present context. Therefore, we shall consider the contribution to the total amplitude coming 
from these stationary points as sub-dominant. This argument is not very rigorous, because in 
this case the spins are not just a scale, but are actually summed over (and doubts could even 
arise about the behaviour of the sum as the cut-off is sent to infinity*''). However, if one think 
to perform the sum first over the spins at fixed scale and then over the different scales,"*^ one 
realizes that the first of these sums should be suppressed at high scales because of the usual 
reason, i.e. because of the absence of stationary points of the action. Nevertheless, the same 
conclusion is supported by the parallelism with BF-theory, where the divergent sector can be 
independently calculated and therefore it is possible to verify a posteriori the absence of the 
contributions stemming from the terms 2 and 3 (see Sect. 3). 

When Sstat = + Hab,a<b ^jabiinfai,, it means that fjab = riab + [i^r]"^, i.e. 

Gab = [±]abGab^SU[2). (9.7) 

''^Recall that the action S appears in the physically meaningful formula only as expS. Therefore, it is defined only 
modulo 2in. 

''^See for example [24] for the discussion of a similar issue. 

''^This would correspond to the change from Cartesian to spherical coordinates in the continuous version of the 
sum over the spins. 
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Here the fact was used that e™°'^ = — 1. Remark that following the discussion at the end of 
Sect. 7.1, the choices of the signs in the previous equations are not independent at every face. 
Rewriting the previous equation in terms of the edge holonomies Ha := gag~^ gives 

Ha = [±]abHb. (9.8) 

With Eq. 9.7 written in this form, it is easy to understand the meaning of the plus or minus 
signs. Indeed, renaming Hi = H, the previous equation just reads 

Ha=±aH, (9.9) 

where it is now emphasized (by the absence of the square brackets) that the ± signs are now 
independent and actually characterize the stationary point. Notice that the holonomy around 
the face [ab) is given by Hab '■= Ha^Hb, and Eq. 9.8 imposes it to be 

Hab=±a±bl- (9.10) 

The solution Hab = I is exactly the solution one would have expected from this sector, which 
is equivalent to the Ponzano Regge model, i.e. three dimensional SC/(2)-BF theory. Neverthe- 
less, with some surprise, it is found that the stationary points treat equivalently ±lsu{2)< and 
are hence only sensitive to the SO (3) structure of the model, viz. to its vectorial representation. 



In the rest of this section we shall calculate the contribution of this family of stationary 
points to the total amplitude. In a first time we shall ignore the presence of possible [m] terms 
in the expressions for the rjab s, which shall later be shovm to be completely uninfluential. 

Equation 9.7 (with the plus sign), implies that 

ga = Kgl and ga = Kgl, (9.11) 

where K,K are arbitrary elements of SL(2,C), and the g° e 5^(2) are functions of the jab'^ 
only. In this notation H = KK~^. Therefore, inserting this result in Eq. 5.2, and following the 
discussion at the end of Sect. 5. 1, we evaluate the partial amplitudes related to this sector to be 

^ji ~A6'^+i2+i2-|-"k(5^sf„\) AKAkW drha {na\Kgl\ma)Y,.jAma\{glr'K-'\na)Y,j. 

JSL{2,Cf^ a JS2 

^^6^+24-irank(<>^Sf,^J d^d^]^ dm„ (n„ |^| m,)Y, j„ (m„ |^| )y, j„ , (9.12) 

where we used the resolution of the identity in terms of the coherent states and the fact that 
the SU{2) group elements commute with the EPRL-FK Yj--map. Notice that the term 6/i in the 
exponent of A comes from the face weights, the first 12 comes from the measures d^- ^rhab, and 
the second 12 comes from the factors i2jab + 1)^ normalizing the measures ^ab^ab- 

This equation shows that at leading order, in the Euclidean sector, the dependence on the 
internal face spins is confined to a scale factor. This fact allows to estimate the divergence de- 
gree of the full amplitude and to explicitly know its dependence on the boundary data. Before 
writing the result, it is important to estimate ran^S^Sfj^j): 

mnk{5^Sfj = 8 x dim(SI(2,C))+ 12 x dim(S2) + 12 x dim(CP^)+ 

eq 

- 4 X dim(SC/(2)) - 2 x dim(SL(2, C)) 
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= 72, (9.13) 

where we calculated as in the BF case of Sect. 3.2.4 the dimension of the space of integration 
minus that of the gauge orbits. Once more, this is only an upper bound, which sets a lower 
bound on the total degree of divergence of this sector. We shall ignore this issue, assuming 
that the bound is saturated as it is in the S[/(2)-BF case. 

Finally, taking into account the six sums over the internal spins, one finds 

pi^Ai?^ ^6(^-1) dKdKYll dma{na\K\rna}Y,jArna\K\na}Yr.ja- (9-14) 

JSL{2,Cf'^ a JS-z 

The possible asymmetry between the solutions of the stationary point equations at the two 
vertices coming from the presence of the [in] 's, does not alter this result. Indeed, at the light of 
the previous discussion, it is clear that one can bring back up the presence of those factors to 
the change in sign of some of the ga 's (or, equivalently, ga 's) 

ga^-ga. (9.15) 
Consequently, the action at the stationary point acquires a term 

Sftat--Sf-,\ + 2i7r^7,fa, (9.16) 

and the amplitudes is modified by an additional factor 

{na\ga\ma)Y,.j, = {na\'irha^Maf^''^{-lf^''{na\ga\ma)Y,,j, ■ (9.17) 

Therefore, the total effect of the change of sign of a ga on the partial amplitude is 

^'^--i-'^f^'-'''"'^',,' (9.18) 

where the convention was used that jaa = ja- However, at each edge of a non-null EPRL-FK 
spin foam graph, the sum of spins must be an integer, since any (formal)*^ integration of the 
type 



dgg>Q^Y^\mi)j, (9.19) 

SL(2,C) i 

subsumes an analogous (well-defined) integration over 517(2): 

dhh>(^\mi)j, (9.20) 

SU{2) i 

which gives as a result an SU{2) intertwiner, and therefore can be non-zero only if 'Y^iii ^ N. 
This fact, combined with Eq. 9.18, concludes our proof. 



''^The reader can check that this interpretation of the origin of the terms [inlab appearing in Eq. 7.22 is coherent 
with the analysis made at the very end of Sect. 7.1. Briefly, if 1712 >-> 1712 + In, this must be traced back to either 
(Si'Si) ^ i~E\'S2) or (gi, g2) "-^ (gi> — §2)- Depending on the possibility that ga and g4 might also have been mapped 
to — gs and — g4 respectively, one finds the different cases listed at the end of Sect. 7.1. 

''^See e.g. [31] and references therein for a discussion about the SL(2,C) intertwiner normalization and SL(2,C)- 
invariance. 



42 SELF-ENERGY OF THE LORENTZIAN EPRL-FK SPIN FOAM MODEL OF QUANTUM GRAVITY 



Therefore, we conclude that the contribution from the Euclidean sector of the sum of the 
internal spins is dominated by a term of order of the form Eq. 9.14. This formula can be 

further simplified by the use of the Dupuis-Livine map Kj- (see Appendix B). Indeed, recall: 

K,{h, g) := I dk ^(27 + lfx^{hk)x^JHkg] (9.21) 

JsUi2) j 



from which it follows 



dg K,ih, g)=\ dg\ dk Jj^j + lfxnhk)x''^Hkg) 

JSL(2,C] JSU{2) j 

dk + ifx'ihk) I d{kg) x'^Hkg] 

j JSL{2,C] 



SL{2,C) 



SUi2) 

= , (9.22) 

where we supposed we could exchange the SL[2, C) and SU[2) integrals as well as the sum over 
the spins, and where the translation invariance of the SL[2, C) Haar measure was used. The 
result follows from the the fact 



dgz"'^(g) = 0. (9.23) 

Therefore, applying the previous formula to the Euclidean dominant contribution Eq. 9.14, 
one gets: 



dgdg]^ dfha {na\g\ma)Y,.jAma\g\na)Y,.ja = 
' a JS9 



SL{2,Cf^ a JS2 

dgdg d/jd/jK^(/z,g)K^[li,g)]^ dma{na\h\ma)jjma\h\ha)j^ 

SL{2,Cf^ JSU{2) a Jsa 

= 0. (9.24) 

9.3 Lorentzian Stationary Points 

The analysis of the contribution of the Lorentzian stationary points to the full amplitude 
closely parallels that given for the Euclidean stationary points. The only relevant differences^ 
comes from the fact that in this case the Gah'^ are genuine elements of SL(2,C) (i.e. ^ SU{2)). 
This means that a crucial step is missing in order to find an equation of the type of Eq. 9.14, i.e. 
the fact that the g^'s do not commute any more with the EPRL-FK Y,.- map. Therefore, the best 
that can be done is to write 

Z ^'^''"'M dKd£]^ I dm,(na|^:g°|m,)Y,j„(m,Kg^)-i^-i|n.)Y,j„, 

{iab<A} JsL{2,Cf^ a JSz 

(9.25) 

where the g^ e SL[2, C) are highly non-trivial functions of the internal face spins {jab}- 



^"Another possible difference could arise from the equations of motions for the "Lorentzian Regge action", which 
may not be the same as for the Euclidean one. However, the expression of the dihedral angles in terms of the spins 
in the two cases is the same, up to an analytical continuation, so the equations of motions read the same in the two 
cases. 
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1 The Total Amplitude 

At the light of the previous calculations, we can conclude that the amplitude of the melon 
graph in the EPRL-FK model, at fixed external coherent states {\na} j^,\na}ja}' ™th arbitrary 
spins [ja], is dominated by a divergent contribution 

^1-^'^""'^ dJCd^]^ dma{na\K\ma)Y„ja{ma\K\na)Y,.u + 

JSL{2,Cf'^ a J $2 

+ ^ A^^/^-^M dKdKYl \ dma{na\Kgl\ma}Y,jAma\iglr'K~'\na}Y,^^^ 

{jab<A} JSLi2,Cf' a JS2 

+ o[A^^t'-'^^). (10.1) 



1 1 About the Approximations 

In this section we briefly review the hypothesis and approximations needed to obtain the 
result presented in the previous section (Eq. 13.6). 

First, we started with the hypothesis that the relevant regime to study the convergence 
properties of the melon graph full amplitude is that in which all the internal spins scale uni- 
formly to infinity. Even if we showed that this hypothesis (when combined to other claims we 
are going to discuss in a moment) leads to the correct result in the context of S[/(2)-BF theory, 
it will continue to be an hypothesis on which all the calculation crucially relies. 

Then, in order to use this hypothesis, we used the stationary phase approximation for the 
integrals. Rigorously speaking this is a viable approximation only if the integration domain 
is compact; therefore in the context of the EPRL-FK model, where multiple integrals are per- 
formed over the non-compact group SL(2,C), a further hypothesis is needed, i.e. some kind 
of well-behavedness of the integrand functions which suppresses their contribution for large 
group elements.^i 

Furthermore, (at least in the context of the non-degenerate sector) we neglected the pos- 
sibility that the Hessian of the action could be degenerate at the stationary points (obviously, 
except for those directions related to the symmetries of the action); to support this hypothesis, 
we showed that it is consistent with our analysis of the self-energy of SU[2)-B¥ theory. 

Finally, we claimed that the terms associated to two tetrahedra with the same parity are 
suppressed with respect to those appearing in Eq. 13.6, because of the strong interferences 
they undergo once the sum over the internal spins is performed; this hypothesis also parallels 
one in the S[/(2)-BF theory, which, in this case, is compatible with to the correct result. 



Some preliminary numerical results seem to show that this is indeed the case. 
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Figure 8: A graphical representation of the "spike" interpretation in four space dimensions. 



12 The Spike Interpretation 

Up to now, we analysed the diverging part of the amplitude in terms of a three dimensional 
"classical" geometry, as our (closure) equations were indicating. In the standard interpreta- 
tion of the divergences, however, they are associated to four dimensional spikes: keeping the 
boundary of the graph fixed (in our cases a topological three-sphere), an additional point is 
added at the "centre" of the triangulation dual to the boundary and is let roam the four dimen- 
sional space. To this roaming is then associated a possible divergence (in the case of BF theory, 
all such configurations are gauge equivalent, and therefore integrating upon them gives a di- 
vergent factor associated to the gauge orbit 

In order to recover this point of view - at least to some extent - we shall consider the "full" 
closure equation at each internal edge: 

eaajama+ ^ eabjabmab = 0, (12.1) 

where eaa = ±1 according to the choice of orientation of the external face a. This equation 
holds once an SU[2) integration is performed at each edge (such an integration is subsumed 
by the SL(2,C) integrations: jsi^^2C]^^"^' ^^^^ equations holds true if and only if the vectors 
{iabh. h^a arid la '■— ^aaja^a <ire the Oriented normals to the faces of an Euclidean tetrahe- 
dron, with length proportional to the face area. 

If the internal spins [jab] are now uniformly rescaled by a factor A ^ oo, while the external 
spins [ja] are kept fixed, the tetrahedra associated to the full closure equations become more 
and more "spiky". Moreover, the normals to the rescaled faces tend to lie on one single plane 
and fulfill Eq. 6.5, instead of Eq. 12.1, to order 0(A~i). Therefore the "reduced" closure equa- 
tion is associated to a normal section of the previous spiky tetrahedron. In this sense, we can 
think we trade a four geometry for a three geometry in our dissertation. 

However, this in not completely correct. Indeed, in order to have a fully consistent four 
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geometry, parallel transport equations must be imposed along the external faces as well; but 
parallel transport equations are an approximate consequence of the largeness of the spins as- 
sociated to a given face. Therefore, nothing of this kind holds for general external faces. 

Nevertheless, if we force ourselves to think in terms of spikes, the natural interpretation of 
the origin of the divergences we found, is the one given in [24], where it was associated to the 
sum over orientation implicit in the EPRL-FK model (as well in most spin foam models). We 
shall come back to this point in the conclusions. 

13 Conclusions 

In this paper we calculated the dominant contributions to the Lorentzian EPRL-FK melon 
graph, i.e. to its self-energy, at fixed boundary states. This is formally done with the insertion 
of a rigid cut-off A on the SU[2) representations characterizing the internal faces of the graph. 
Then, in order to study the convergence properties of the sum, the relevant large-spin limit of 
the partial amplitudes has been taken. This allowed us - under some hypothesis - to use the 
stationary phase approximation for the integrals. Hence, the equations defining the station- 
ary points have been solved explicitly, and the solution have been interpreted in geometrical 
terms. The dominant contributions, for those sets of spins which allow their existence, are 
found to be given^- by two equal Euclidean, or Lorentzian (depending on values of the spins 
one is considering), tetrahedra of opposite parity. 

The result about the amplitude is then twofold. On one side we calculated the self-energy 
scaling in A, as a function of the spin-foam face-weight [lij). On the other side, we worked 
out the structure of the dependence on the boundary states of the dominant part of the graph 
amplitude. 

The scaling of the EPRL-FK self-energy is found to reflect the Euclidean results of [9, 8].^^ 
In particular the dominant contribution diverges as 

Wj^ 0[A^^''-^'>) (13.1) 

in the spin cut-off A. Here, [i is defined by the asymptotic behaviour of the spin-foam face 
weight: /i(A) ~ A^^ as A ^ oo. In particular, this means that if the face weight is the one imposed 
by the invariance of the spin-foam amplitude under face splitting {^[j] = 2j + 1), then the 
dominant (and consequently only) divergence present in the melon graph is logarithmic in 
the cut-off: 

fjiU) = 2j + l ^ Wj^^OilnA). (13.2) 

At this point one could try and interpret the cut-off model as a qualitatively viable and more 
treatable version of the full ^-deformed model. In the ^-deformed model a cut-off on the S[/(2) 

^^With the caveafs thoroughly discussed in Sect. 11. 

^^This is done in a different setting from that used in [9, 8] (though the one of [8] is similar). Nonetheless, our 
calculation is straightforwardly generalizable to the Euclidean case and, at the light of this fact, it is clear that the two 
results must coincide. Indeed, our technique essentially calculates the dimension of the integration space used to 
calculate the amplitude minus the dimension of the symmetries of the amplitude itself. Moreover the only relevant 
variables for the calculation turn out to be the group variables (this is not explicitly shown in the text, but the reader 
can rapidly demonstrate it by showing that the contributions coming from the "auxiliary" variables [rha] and [Zat] 
cancel thanks to their measure scalings). Finally, since dimS0(4) — dimSL(2,C), the equality between the Euclidean 
and Lorentzian results follows. 
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representations arises naturally, and clues exist which point to the relation of this cut-off with 
the cosmological constant Acc- In particular, A ~ Acc/lp ^ lO^^o, Ip being the Planck length. 
Then the logarithm appearing in Eq. 13.2 would not be very large at all! Indeed: 

ln/l~ln(10^20)^280, (13.3) 

which is a small number when compared to those one would expect when considering repre- 
sentations of the order of A itself. 

Moreover, the previous estimation seems to be meaningful only for radiative corrections of 
phenomena taking place at the Planck scale: if one takes into account the fact that there exist 
a scale at which our analysis in terms of the stationary point approximation starts to be viable, 
and that this scale is set by {some large number times) the boundary spins,^* then Eq. 13.2 
should be replaced with something more similar to 



M^^.~Olln^l. (13.4) 



V Jext 

Hence, for physics at the galactic scale jext ~ (10^1y)^/Z| 10^^". This implies 

ln^~ln(10^°)?«20, (13.5) 
Jext 

which is essentially of order 1!^^ Therefore, this suggestive setting leads us to speculate that 
in a theory with cosmological constant, the supposedly very large divergences associated to 
melon graphs could be, in the end, not large at all (but not small either!). In general, the result 
would clearly be scale dependent, and the larger the (length) scale, the smaller the corrections. 
Remark, however, that the framework we used for the calculations breaks down at the cosmo- 
logical scale, since it is strongly based on the existence of a neat scale hierarchy: jext <■ A. 

For what concerns the other result, i.e. the dependence of the dominating part of the am- 
plitude on the boundary states, it reflects the very nature of the EPRL-FK model. Indeed, its 
form strongly relies on the EPRL-FK map. We report here the final formula: 

^i-^'^'^'^M d^:d£]^ dma{na\K\ma}Y,.jAma\K\na}YrJa + 

JSL[2,Cf^ a JS2 

+ d^d^n I dm,(n„|^g«|m,)Y,j„(m„|(g°)-i£-i|n„)Y,j„^ 

{iab<A] JsL(2,Cf2 a JS2 

+ o{A^^i'-^'>). (13.6) 

The second term in this formula is a new feature of the EPRL-FK model with respect to BF- 
theory. It stems from the presence of a Lorentzian sub-sector of the solutions to the stationary 
point equations. In this sector the group elements g'^ are genuine SL(2,C) elements and do 
not commute with the EPRL-FK Y,, map. The explicit evaluation of this term is obstructed by 
the implicit (and very intricate) dependence of the on the internal spins {7a/,}.0ne can try 
and guess that because of this dependence, the sum over the internal spins gets suppressed. 
This would give a result quite intriguing for its simplicity: 

^z-^'^""" d^d^f^ dma{na\K\ma)Y,,jAma\K\na)Y,.u + o{A''^''-'^). (13.7) 

JSLiZ.Cf^ a J $2 



^Otherwise the closure equation in the form of Eq. 6.5 would not be accurate. 



^^At our scale jext ~ Im^/ll ^ 10''^, hence In -A- ~ln(10^^)f^ 115, which is of the same order as Eq. 13.3. 
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The structure of this formula reminds the EPRL-FK bare propagator, which coincides with that 
of Sf/(2) BF-theory.56 Formally, 

P:= dh0h, {na\ma}u^nja,na,na)=\ dhY\{na\h\na}u. (13.8) 

JSU{2) a JsUi2) a 

To stress the similarity of Eq. 13.7 with the propagator, and to write it in a more compact form, 
define the (formaP^) operator T^: 

T^:= d^(g)Y/^^, {na\Tj\na)j„^Tj{ia,na,na)=\ AKY[{na\K\na\.ja- 

JSL(2,0 a JSL{2,C) a 

(13.9) 

Remark that T,, is morally nothing else than "K^ > P" (see Appendix B). Notice, also, that 
= PTj- = TjF, since the SL(2, C) integration subsumes an SU{2) integration. In terms of this 
new operator, the EPRL-FK self energy at its dominant order (Eq. 13.7) reads simply 

Wj^[ia,na,na)^A''^^'-'^\na\'T\\na)u ■ (13-10) 

Remark that is not a projector (on the contrary P is), therefore ^ Tj. This is because 
of the presence of the EPRL-FK Yj map and is therefore intrinsic to the nature of the model. 
For this very same reason, T,. is an operator which appears in many practical applications of 
the EPRL-FK model (e.g. in Spin Foam Cosmology [41]) and therefore deserves a study of its 
general properties. For example, its semi-classical behaviour at large spins {ja} was recently 
unravelled in [42] . 



Finally, remark also that the most diverging part of the (non-degenerate sector of the) melon 
graph amplitude is not proportional to the bare propagator of the model. The consequences 
of this fact on the renormalization group of the EPRL-FK model and on the physical implica- 
tion of the model itself rest to be studied. Nonetheless it seems clear that the dynamic of the 
renormalized model would result modified with respect to the tree level one. Whether this 
modification goes in the direction of invalidating the results about the EPRL-FK semiclassical 
limit, or rather in the direction of enriching them (e.g. by modifying the semiclassical limit 
to an approximation of General Relativity richer than Regge calculus), is still an open ques- 
tion. For the moment, we just remark that a possible mechanism for this is the emergency of 
non-trivial gluing conditions due to the dependence of the renormalized propagator on the 
spins {ja\- Indeed, something very similar happens in the S[/(2)-BF case (see [36]), where the 
result of the melon renormalization procedure is the introduction of a group Laplacian in the 
renormalized GET Lagrangian. 
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Problems could arise from the non-compactness of SL{2, C). 
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A Euclidean or Lorentzian? 

In this section we shall give a more straightforward characterization in terms of the spins of 
Euclidean or Lorentzian nature of the stationary points. Indeed, in Sect. 7.1 a formula is given 
which allows for a distinction of the Euclidean from the Lorentzian case, in terms of the spins 
(Eq. 7.16). However, now that we have a correspondence between the spinorial data at the 
stationary points and geometrical tetrahedra, we can analyse the same issue from a different 
perspective. 

In order for the £ab's to be solutions of the closure equations, the spins {jab}b,b^a must 
satisfy triangular inequalities. We shall suppose these inequalities hold strictly, i.e. we shall 
not consider the case of degenerate triangles. One can ask if further conditions are needed to 
define a tetrahedron in M^'^. To answer this question, we go through the following geomet- 
rical construction. It was shown that all the faces of the considered tetrahedra are spacelike. 
Therefore, it is always possible to boost the whole tetrahedron in such a way that face 4 lies on 
the xz -plane. Further rotating the tetrahdron around the y axis, allows to align side 13 to the 
X-axis. Finally, translating the tetrahedron one can make coinced vertex 1 with the origin of 
the Cartesian system (see Fig. 9) . In this configuration, the equations for the tetrahedron edge 
lengths read 





— X3 


% 


= X2 + ^2 


J23 


= {X3-X2f + Zl 




= xl + zl + yl-tl 


% 


= (X4 -X3f+zl + y| - f| 


. -^24 


= (X4 - xzf + [Zi - zzf + yl - tl 



(A.1) 



where {xi,yi,Zi, ti) are the four coordinates of the vertex / = 1, . . .,4. The first three equations 
concern the triangle 1 and involve the resolution of a quadratic equation, whose solution is 
guaranteed to exist by the triangular inequalities. Subtracting equation four from five gives a 
linear equation for X4, and so does the subtraction of equation five from six for Z4. Thanks to 
their linearity, these equations are readily solved in terms of the j 's and of the solution to the 
first three equations. Equation four becomes then an equation for [z^ - f|): 

yl-tl=il^-xl-zl^Q, {A.2) 

where the right hand side of this equation can be either positive or negative. Provided that 
the z f -plane is appropriately boosted, in the first case the tetrahedron lies in the hypersurface 
orthogonal to and is said to be Euclidean; while in the latter case, it lies in the hypersurface 
orthogonal to and is said to be Lorentzian. 

Remark that during the previous procedure, two quadratic equations had to be solved. This 
implies that two sign choices had to be made. The first of these choices corresponds to the 
choice of orientation of the tetrahedron basis (the face lying in the xz -plane). The second one 
correspond to the sign of 3/4 or in the Euclidean or Lorentzian cases, respectively. Therefore 
any set of six edge lengths {jab} satisfying triangular inequalities for any subset {jab}b,b^a, 
defines up to global rotations and boosts, and up to space and time inversion, a tetrahedron in 
M}'^. The tetrahedron can be either Euclidean or Lorentzian. In the first case the edge lengths 
satisfy also the condition of positivity of the so-called Caley-Menger determinant [43, 24] . 
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Figure 9: The six edge lengths {ji2,ju,j23,ju,j2A,j34] closing into a tetrahedron in M*, with the axis 
chosen as in the text. 



B The Dupuis-Livine Map and Face Reversal 



In this appendix the following equation is shown to hold, from which it is straightforward 
to deduce the face reversal symmetries of the EPRL-FK model: 



Tr; 



= Tr,- 



(B.l) 



In order to to this, it is useful to introduce the Dupuis-Livine map [44] K^(h, g): 



^^[51(2,0] 



(B.2) 



whose kernel, with litde abuse of notation, is defined by 

K,{h,g)~ I dfc^(2j + l)2j^(^fc)j"J(A:g). 

JsUi2] j 



(B.3) 



Here and /P *^ are characters of SU{2) and SL(2,C) in their j and [p, k) representation re- 
spectively. 

From this appendix's point of view, the crucial property of Dupuis-Livine map is that 

K^{h,g) = K^{h-\g-'). (B.4) 

To show this, beside the fact that jgjj^^^)'^^ ~ hu{2)'^^~^' enough to recall that xK^) = 
X-'ih'^) and x'^'^'ig) = The first of these identities follows from h = h~^^ and the 

fact that SU[2) characters are real,^*^ while the second one can be shown by mean of the map 



ssExplicitely, using the fact that V/i e SU[2) 3ue SU[2), 6 e [0, 2n] such that h^u e^'^"- u-^: 



m=^j sin [2^] 



50 SELF-ENERGY OF THE LORENTZIAN EPRL-FK SPIN FOAM MODEL OF QUANTUM GRAVITY 



[38, 31, 11] / : ^P'*^ ^P-*^. Indeed, it has the following properties: 

i/g^>g^%,k = {Jip,ip%.k and iip,ip%,k = {Jip',f4>)p,k Vj/-,j/;'e=^P'*^VgeSi(2,C) 

(B.5) 

where (•, Op./t stands for the Hermitian scalar product in ^p-*^. From these, it follows that 

if^>gf^)p.k = [fg-'ip.ip)p,k = [ip,g-'ip)p.k Vj/'e^^P'^ (B.6) 

Therefore, since if is an orthonormal basis then {Jipi}i is an ortho normal basis as well, 
by summing the previous equation over such a basis, one immediately finds the sought con- 
clusion: x^Hg) = x'''^ig~^)- 



Using Schur's orthogonality relations, Eq. B.3 can be rewritten as 



(B.7) 



where no sum is taken over the indices m and n, and where D-i[h) and DP'^[g) are the ma- 
trices of the elements h e SU{2) and g e SL(2,C) in representation j and {p,k), respectively. 
Therefore, using Peter- Weyl theorem, one finds 



(B.8) 



Formally being Drhi{g)>^^ = {j ,m\Yj^ Drhi{g)Yj\i , n), this can be thought as 

<?(g):=(K^>0Xg)^(^(Y/gY^). (B.9) 

Using this notation, it is immediate to see that a face (partial) amplitude in the EPRL model 
is the transformed through of an SU(2)-BV one: 



Tr; 



e \ e JSU[2) \ e J 



(B.IO) 



Now, the partial amplitude of the same face of the previous equation, taken with opposite 
orientation is 



Tr,- 



= Y\\ dHe K,(//„ G;1) x'' Yl^e 

e JsUi2) V e J 

= dH,KriH;\Ge)xjiYlH;' 

e JSU{2) V e 



= Tr,- 



(B.ll) 



where Eq. B.4 was used, as well as the fact that J^^j^g) ~ hu{2] ^ ^'^^ X^W = J ^ (^ ^) for 
heSU{2]. 
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C The Degenerate Sector 

In Appendix C of [8], it was siiown that the geometrically degenerate sector of the melon 
graph in the Euclidean EPRL-FK model contributed to the graph amplitude more than the non- 
degenerate one. This result was achieved via explicit estimation of the asymptotics of the 15j 
and 9 j symbols building the Euclidean version of the EPRL-FK vertex amplitude. Moreover, 
this was done by fixing the external spins (the [ja } in our notation) exactly to zero. 

Here, we try to translate their result to our language, which better suits the Lorentzian ver- 
sion of the model. This shall also clarify the geometrical origin of the dominance of this sec- 
tor in the EPRL-FK model, while it is 5U;b-dominant in the Sf/(2)-BF theory (as was indirectly 
shown in Sect. 3). Nevertheless, we shall not go into the details of explicitly evaluating the con- 
tribution of this sector. 

We start by classifying the degenerate sectors of the melon graph. In particular, there exists 
three subsectors of degenerate configurations, depending on how many triangular inequalities 
are saturated: 

1. the triangular inequalities are strictly satisfied, but the Caley-Menger determinant of the 
spins is null; 

2. at one (and only one) edge the triangular inequalities are saturated; 

3. at more than one edge the triangular inequalities are saturated . 

This sectors correspond to geometrical configurations in which the tetrahedron is flattened to 
lie on a plane in three different ways: 

1. all the four vertices lie on a plane, but no numerical coincidence exists among the vertex 
positions (the faces are triangles with non-zero area); 

2. a face is degenerate, therefore three vertices lie on a line, while the fourth one can lie 
anywhere in space, but on the line spanned by the three other vertices. 

3. the four vertices are aligned. This sector shall be dubbed maximally degenerate. 

For the previous coincidences to be possible, a certain number of constraints must be satis- 
fied by the [jab] in such a way that the six independent sums on the j 's are effectively reduced 
by 1,2 and 3 sums, respectively. At the same time, the stationary points acquire more sym- 
metries; this fact tends to enhance the divergence related to this sector. The balance between 
these two phenomena shall determine the relevance of the degenerate sector over the non- 
degenerate one. 

For simplicity in this section we will only analyse the maximally degenerate case, i.e. case 
3. here above. Moreover, we will also avoid to go into the algebra of the stationary point equa- 
tions, even if this would be almost straightforward in the maximally degenerate case. There- 
fore, the presented arguments will be essentially of a geometrical nature. 

For a better understanding, we start from the simple case of SU{2)-B¥ theory. At each edge, 
the degenerate configuration is represented by three points aligned on a segment. It is intu- 
itively clear that a new symmetry must associated to this configuration. Indeed, this is the 
case, and the symmetry is a U{1) symmetry around the direction of the edge itself. Therefore, 
the four parallel transports within a tetrahedron ha (and ha) are all determined up to one of 
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these [/(I), except for the fact that a constraint among these symmetries exists. One can re- 
alise that this is the case by observing that in the spinor representation aU[l) rotation along 
the spinor axis introduces a phase, and the four phases so introduced within a tetrahedron 
must compensate each other. Finally, four of these symmetries can just be reabsorbed in the 
SU{2) rotational freedom we have at each edge (i.e. can be reabsorbed in the definition of the 
edge reference frames). As a result the total rank of the Hessian is (at most) reduced of 

rank(52s»Pg) - mnki5^Sf^^_^^^] = -(8 - 2 - 4) x dim([/(l)) = -2. (C.l) 

Therefore, taking into account that there are 3 less free spins on which to sum, one can esti- 
mate the degree of divergence of the degenerate sector in the S[/(2)-BF theory: 

Wf, ~yl-3xyl^H'^^^''''s)-''''"'^(^''n-^^^ , ~A-2x</ , , (C.2) 

^(,aeg ^/i', non-deg .-«, non-deg' ^ ^ 

which means that the degenerate sector is subdominant with respect to the non-degenerate 
one. 

Unluckily, things might different for the EPRL-FK model. The reason is that in this model 
the vertex boundary states are the same as in SU[2)-BF theory, however the parallel transports 
double the copies of SU[2) present in the bulk.^^ As a consequence, the symmetries on the 
degenerate boundary states induce a larger freedom in the bulk. For clarity we shall treat the 
Euclidean and Lorentzian EPRL-FK models independently, even if they lead to the same result. 

In the Euclidean model, each Spin[4) holonomy is made of a self-dual and an antiself- 
dual part g = ig'^,g~), each in SU[2), and whose combined action on the boundary states 
is dictated by the fusion coefficients (see [3]). Therefore the counting we did in the SU[2)- 
BF theory applies almost independently on the self-dual and antiself-dual, the only exception 
being the fact that the freedom to redefine the edge reference frame is not doubled (the edge 
boundary state being the same as in the previous case!). Therefore, 

mnk[5^S^f) - rank(52sE-i^^g) = -[2 x (8 - 2) - 4] x dim([7(l)) = -8. (C.3) 

Hence, 

W^]""} ~A-3xA^['''"''^^'^^eg')-'-an«52S^-l,^g)] ^ E^^^^ ^^1 ^Eucl 

^,deg non-deg .41 , non-aeg ^ ' 

This result is in complete accord with the calculation performed in Appendix C of [8] , with the 
only caveat that they are using a different face weight''^ (essentially the BF theory one, which 
has/i = 2). 

In the Lorentzian case, being the parallel transports in SL(2, C), thanks to the symmetries of 
the degenerate boundary state, they are defined only up to f/(l)c element, which has dimen- 
sion 2!°^ On the other side, only 4 real can be reabsorbed in the redefinition of the edge 



^^This is exactly true in the Euclidean EPRL-FK model, and morally true in the Lorentzian one, where the 
holonomies are in SL(2, C), whose algebra, in the fundamental representation, is essentially a complex, instead of 
real, St/(2) algebra. 

s^Remark that the GET construction is less flexible in the choice of the edge and face weights, than that of Spin 
Eoams. 

^^The physical interpretation of this is related to the possibility of boosting the faces along their degeneracy direc- 
tion. 
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reference frames. Therefore, 

rank(52s^^g^') - rankCS^S^^;;!^^^) = -(8 - 2) x dim([7(l)c) - 4 x dim([/(l)) = -8. (C.5) 
And finally, 

WC#,deg ~ X Ai [rank(5^S,eg)-rank(5^S„o„-<ieg)] x l^^.^on-deg ~ vl^ X W^,non-deg. (C.6) 

In conclusion, this estimation hints to the fact that the geometrically degenerate sector 
contributes more than the non-degenerate one to the EPRL-FK self-energy. Nevertheless, one 
has to be cautious. Indeed, the degenerate sector happen to be close to the boundary of the 
region of validity of some of our assumptions and the treatment we gave to the problem could 
break down and must be carefully checked. On the other side, this sector might be simple 
enough to allow an explicit analysis of the Hessian of the action. 
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